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1.  Summary  of  Progress 


Work  has  proceeded  according  to  schedule  on  all  research  topics. 

(i)  Uniquely  well  defined  clusters  of  Gold  are  now  available  for 
optical  measurements. 

(ii)  A  comprehensive  experimental  investigation  has  been  carried 
out  to  determine  the  non-linearity  of  these  samples  on  a 
femtosecond  timescale.  These  results  are  presented  below  and 
prove  that: 

(a)  the  nonlinear  refractive  index  rises  and  relaxes  on  a 
sub-picosecond  timescale. 

(b)  the  nonlinearity  is  predominant  positive  imaginary  originating 
from  a  hot  electron  contribution. 

(c)  the  nonlinearity  is  independent  of  the  size  of  the  particle. 

2.  Staffing 

No  major  changes  in  staffing  have  taken  place. 


CHAPTER  1 


INTRODUCTION 


The  interaction  of  light  with  matter  is  usually  characterised  by  several 
phenomena,  such  as  light  absorption,  refraction,  scattering  etc.  All  of  these  are 
regarded  as  linear  optical  properties  of  the  material,  dependent  on  wavelength  but 
independent  of  the  intensity  of  the  light.  However,  if  the  illumination  is  made 
sufficiently  intense,  the  optical  properties  begin  to  depend  on  the  intensity  and  other 
characteristics  of  the  light.  The  light  waves  may  then  interact  with  each  other  as  well  as 
with  the  medium.  The  study  of  such  interactions  is  the  field  of  nonlinear  optics.  It  is  a 
field  that  has  experienced  a  rapid  growth  of  interest.  There  are  three  main  motivations 
for  this.  First,  there  is  the  possibility  of  exploiting  the  nonlinear  behaviour  in  various 
devices.  The  most  important  of  these  are  frequency  converters,  in  which  laser  radiation 
at  one  frequency  is  converted  into  coherent  radiation  at  a  new  frequency,  by,  for 
example,  harmonic  generation,  sum-frequency  generation  ele.H)  Because  the  converted 
radiation  may  be  at  a  frequency  that  is  not  directly  available  from  a  laser  source,  these 
frequency  conversion  techniques  provide  an  important  means  of  extending  the  spectral 
range  covered  by  coherent  sources.  A  second  reason  for  studying  nonlinear  optical 
processes  is  that  they  set  a  limit  to  the  light  flux  that  can  be  passed  through  a  medium. 
For  example,  two  photon  absorption  can  lead  to  depletion  of  the  incident  light  and  self- 
focusing  leads  to  distortion  of  the  incident  beam  profile.  A  third  interest  in  nonlinear 
optical  effects  lies  in  their  use  as  a  means  of  obtaining  information  about  the 
microscopic  properties  of  the  atoms  of  molecules  that  constitute  the  nonlincr  medium. . 
Two  photon  absorption,  for  example,  can  be  used  to  study  energy  levels  that  are 
inaccessible  by  single  photon  absorption. 

When  an  electric  field  is  incident  on  a  conductor  consisting  of  negatively  and 
positively  charged  particles,  it  gives  rise  to  a  flow  of  current  -  the  positive  charges 
move  in  the  direction  of  the  field  while  the  negative  ones  move  in  the  opposite 
direction.  In  the  case  of  dielectric's,  on  the  other  hand,  the  charged  particles  are  not  free 
to  move  but  arc  bound  together.  The  effect  of  an  electric  field  on  such  materials  is  to 
induce  a  polarisation  that  is  proportional  to  the  applied  field.  However,  as  the  intensity 
of  the  field  is  increased  the  linear  dependence  of  the  two  terms  no  longer  holds.  The 
polarisation  must  now  be  expressed  as  a  power  scries  in  the  field,  as  follows  : 

P  =  en(XwE+xa)L T2+*l3’£3+ )  ,,  n 
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Where.  >  Linear  su.scepiihilitv.  aiul 
Nonlinear  susceptibilities. 

The  nonlinear  .susceptibility  is  a  tensor  object  and  reflects  the  symmetry  of  the 
material.  l;or  a  crystal  with  inversion  symmetry',  i.e.  one  that  is  centrosymmeiric,  the 
second  term  and  higher  even  terms.  etc.,  are  zero.  However  %0)  is  always 

non  zero. 

At  first  crystalline  nonlinear  media  were  seen  to  offer  the  greatest  device 
potential.  The  reason  for  this  lay  in  die  fact  that  %l2)  is  zero  unless  the  medium  lacks  a 
centre  of  symmetry.  Nonlinear  effects  in  centrosymmeiric  media  would  thus  depend  on 
higher-order  and  therefore  presumably  smaller  nonlinear  terms.  Unfortunately, 
however,  few  crystalline  materials  have  proved  capable  of  satisfying  the  list  of 
requirements  for  a  good  nonlinear  material.  A  number  of  significant  advantages  are 
offered  by  liquids  and  gases.  They  can  be  easily  prepared,  with  good  optical  quality 
over  large  dimensions.  Also,  they  are  less  prone  to  suffer  irreversible  damage  at  high 
intensities. 

In  order  for  nonlinear  behaviour  to  manifest  itself,  the  field  incident  on  the 
material  must  be  comparable  to  the  internal  field  E,  which  binds  together  its  electrons 
and  ions;  typically  Ea~  3  x  1010  Vm'1 .  To  obtain  an  optical  field  of  such  a  magnitude, 
an  incident  intensity  of  =  10U  Went'1  is  required.  However,  such  high  intensities  are 
not  in  fact  necessary  for  the  observation  of  many  nonlinear-optical  effects.  One  reason 
is  that,  provided  the  assembly  of  induced  dipoles  oscillates  coherently,  the  field  that 
they  radiate  individually  can,  in  certain  circumstances,  add  together  constructively  to 
produce  a  much  larger  total  intensity.  This  phenomenon  is  termed  "phase  matching". 

The  intensity  required  to  observe  some  nonlinear  processes  can  be  further 
reduced  by  many  orders  of  magnitude  by  choosing  one  or  more  of  the  optical 
frequencies  involved  so  that  they  lie  close  to  a  resonant  frequency  of  the  oscillating 
dipoles;  this  is  termed  "resonance  enhancement".  In  nonlinear  optics  resonant 
enhancement  is  utilised  in  two  ways ;  Firstly,  it  allows  nonlinear  processes  and  devices 
to  operate  effectively  at  lower  power  levels,  thus  increasing  their  range  of  use  and 
efficiency.  Secondly,  resonant  nonlinear  phenomena  provide  the  basis  for  "nonlinear 
spectroscopy";  the  observation  of  these  effect  can  provide  information  about  the 
structure  of  matter  that  is  not  accessible  using  conventional  linear  optical  spectroscopy. 

In  order  to  be  useful  in  an  optoelectronic  device,  a  material  requires  a  large 
nonlinear  susceptibility  along  with  a  last  response  time.  The  aim  of  this  study  is  to 
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determine  what  properties  ol  a  material  come  imo  play  in  the  enhancement  ol  these 
parameters. 


THE  SUSCEPTIBILITY  TENSOR 


By  considering  the  dynamic  behaviour  of  charged  particles  in  a  medium  under 
the  influence  of  an  electric  field,  it  is  possible  to  derive  explicit  formulae  for  the 
susceptibility  tensors  ol'  the  medium.  For  simplicity  it  is  assumed  that  the  medium 
consists  of  an  assembly  of  microscopic  polarizable  units  which  are  independent  and 
non  interacting.  Thus  local  field  factors  are  neglected.  Local  field  corrections  can  be 
introduced  later  without  altering  the  general  formuL 

A  straightforward  derivation,  using  the  density  matrix  approach  leads  to  the 
following  formula  for  the  third  order  susceptibility  tensor  £<3>,  that  depends  on  the 
frequencies  6),,  co: ,  and  <y,  of  the  three  electric  fields  involved 
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Where  fi,  a,  fi,  y  arc  summed  over  the  co-ordinates  x,  y  and  z, 

(Og  =  <y,  +  co2  +  cov  where  <yp  co2,  and  tu,  are  the  frequencies  of  the  applied  fields, 
the  operator  S  implies  intrinsic  permutation  symmetry, 

N  is  the  number  density  of  molecules,  i.e.  the  number  of  molecuies/unit  volume, 
er ^  is  the  (ab)^  element  of  the  dipole  moment  matrix  er* ,  and 
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E  -  E 

— - —  is  the  transition  frequency. 
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This  expression  has  singulai  ines  in  t Ik-  Irequencv  domain  whenever  any  oik-  o! 
the  frequency  denominators  of  the  type  L2  -  o)  approaches  zero.  Thus  the  susceptibility 
may  be  "resonantly  enhanced"  by  choosing  one  or  more  of  the  incident  frequencies,  or 
combinations  of  them,  to  coincide  with  an  optical  transition  frequency  of  the  medium. 
Obviously,  as  is  apparent  from  the  above  formula,  for  the  largest  nonlinear 
susceptibilities  we  require  large  values  for  the  matrix  elements  of  the  electric-dipole 
operator  cr  and  small  frequency  denominators. 

The  use  of  resonance's  to  enhance  a  specific  susceptibility  often  has  the  effect 
of  similarly  enhancing  other  competing  or  undesirable  processes.  For  example,  by 
allowing  an  input  frequency  to  approach  resonance  with  a  single-photon  molecular 
transition,  that  input  field  may  suffer  a  corresponding  increase  in  absorption;  this 
occurs  because  the  first-order  susceptibility  £(l)  is  also  enhanced. 

For  very  close  resonance,  the  mathematical  divergence’s  in  the  above  formula 
are  unphysical;  they  occur  only  because  higher  order  nonlinearilies  have  been 
neglected.  When  excited  very  close  to  resonant  transitions,  the  molecules  undergo  large 
perturbations  -  thus  invalidating  the  small-perturbation  approximations  -  and  the 
transition  frequencies  which  occur  in  the  denominator  terms  become  field-dependent 
themselves.  When  these  strong  field  dependent  perturbations,  or  level  shifts,  are  taken 
into  account,  the  resulting  induced  polarisation  remains  finite.  In  many  cases,  however, 
the  resonant  nonlinearities  are  dominated  by  various  transition-line  broadening 
processes,  perhaps  due  to  interactions  between  the  molecules,  which  also  ensure  that 
the  resonant  susceptibilities  do  not  diverge.  A  particular  resonant  process  can  then  be 
represented  by  a  single  order  of  nonlinearity  derived  using  the  small-perturbation 
analysis,  but  with  the  addition  of  appropriate  damping  terms  ±iT  in  the  frequency 
denominators  of  the  susceptibilities  as  for  example  +  iFhj  -  ioa).  The  damping 

factor  Fha  is  thus  identified  as  a  dephasing  parameter  appropriate  to  weak  collisions, 
and  as  such  represents  a  spectral  linewidth  for  the  transition  of  frequency  Qba.  The 

smallest  values  which  the  damping  factors  T  can  take  arc  determined  by  spontaneous 
emission,  and  in  this  case  rb~'  is  a  natural  lifetime  in  the  absence  of  collisions  and 
other  perturbations. 

Thus  when  damping  terms  arc  inserted  into  the  expression  for  the  general  third- 
order  susceptibility,  the  following  expression  is  obtained  : 
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Clearly,  if  ihe  frequency  which  appears  in  a  particular  denominator  term  is  far 
removed  from  molecular  resonance's,  then  the  corresponding  damping  term  can  be 
safely  neglected. 


A  price  paid  for  the  resonant  enhancement  of  optical  nonlinearity  is  the  slowing- 
down  of  the  speed  of  response,  which  may  he  a  serious  limitation  for  some  applications 
such  as  optical  switching  and  signal  processing. 


The  work  done  for  this  thesis  is  based  on  the  measurement  of  the  resonantly 
enhanced  nonlinear  susceptibility,  of  colloidal  gold  solutions,  with  the  diameters 
of  the  gold  particles  varying  from  50-4(X)A.  The  reason  for  interest  in  these  materials 
lies  in  the  fact  that  the  properties  of  the  inclusions  are  completely  different  from  hulk 
gold  because  of  their  small  dimensions.  This  results  in  two  confinement  effects, 
quantum  confinement  and  dielectric  confinement,  which  both  have  potential  in 
enhancing  the  macroscopic  susceptibility  of  the  solution  and  are  also  in  themselves  very 
interesting  phenomena  to  study.  In  chapter  2  the  materials  are  described  in  more  detail 
and  the  confinement  effects  mentioned  here  arc  elaborated  upon.  Chapters  3, 4, 5  and  6 
describe  different  experimental  techniques  that  were  used  to  study  the  materials  and  in 
particular  to  measure  their  susceptibility  and  response  times.  Chapter  7  gives  the  final 
conclusions  of  this  study  based  on,  both,  the  results  obtained  here  and  findings  in  the 
extensive  literature  on  this  topic. 
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CHAPTER  2 


THH  OPTICAL  PKOPKRTIES  OF  ARTIFICIAL  MATERIALS 


In  the  early  Jays  ul  nonlinear  opues.  the  materials  used  lor  experiments  and 
devices  were  mainly  inorganic  dielectric  crystals,  vapours,  liquids  and  hulk 
semiconductors.  More  recently,  interest  has  been  focused  on  new  artificial  solid  state 
materials  which  offer  higher  nonlinearity,  and  in  particular,  those  that  allow  nonlinear 
optical  devices  to  operate  el  I  lciently  at  relatively  low  power  levels.  Several  materials 
with  large  second-order  nonlinearity  have  been  successfully  fabricated  1  "These  have 
applications  in  devices  such  as  compact  optical -frequency  doublers  and  parametric 
amplifiers  and  oscillators.  However,  materials  with  large  third  -order  nonlinearity  are  ol 
greater  interest  currently,  since  the  nonlinear  refracuve-index  effect  can  be  exploited  lor 
switching,  optical  bistability,  phase  conjugation  and  other  types  of  signal  processing1-'. 
A  lot  of  advances  have  been  made,  of  late,  in  the  synthesis  of  organic  materials, 
inorganic  semiconductors  and  metallic  clusters  and  particles  with  a  large  third  order 
nonlinearity.  The  emphasis  of  this  study  is  on  composite  materials  composed  ol  small 
metallic  particles  embedded  in  a  dielectric  matrix.  A  composite  medium  is  made  of  tat 
least)  two  different  component  media,  but  it  is  not  an  alloy  as  they  are  not  mixed  at  the 
atomic  level. 

The  linear  optical  properties  of  metal  colloids,  consisting  of  metal  particles 
suspended  in  a  dielectric  matrix  have  been  studied  for  a  long  limc/T  4)  The  first  models 
aiming  at  describing  these  properties  were  developed  at  the  turn  of  the  century.0- 
However,  it  was  not  until  the  1980’s  that  these  materials  were  studies  as  nonlinear 
optical  materials.*7- Subsequently  a  lot  of  work  has  gone  into  investigating  them,  ilie 
discovery  that  they  exhibit  a  relatively  iarge  oand-edge  resonant  third-order  nonlinearity 
with  relaxation  times  as  short  as  5ps  has  stimulated  much  of  the  current  investigation  'T 
10,  11,  12,  13)  -file  linear  dimensions  of  the  conducting  particles  in  these  composite 
materials  are  of  the  order  of  a  few  nanometers  and  the  carriers  in  them  arc  confined  in 
all  three  directions.  It  is  therefore  customary  to  refer  to  them  as  quantum  dots  or  zero 
dimensional  structures  *14'.  This  term  is  appropriate  only  when  the  dimensions  of  the 
confinement  arc  smaller  than  the  mean  free  path  of  the  carriers.  Then  the  description  of 
the  optical  behaviour  involves  the  energy  eigenvalues  and  eigenfunctions  resulting  from 
the  carrier  confinement. 
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2.1  MATERIAL  PREPARATION 

The  materials  used  in  this  study  consisted  of 'colloidal  suspensions  of  different 
si/ed  cold  particles  in  deionised  water.  The  diameters  of  the  cold  panicles  ranged  1n>m 
S0-400A.  The  sols  arc  labelled  C(50A).  K(172A).  M(592A).  N(.W5A).  Ot.ttrfA). 
P(52A).  Q(47 A).  R(32f>A).  Y(179A).  Z(164A)  and  SflSdA).  Sols  K.  M.  N.  ().  R.  V 
and  Z  were  prepared  according  to  the  method  of  Frens  *'“0  -  reduction  of  chloroauric 
acid  solution  bv  Sodium  citrate,  which  (when  successful)  produces  possibly  the  most 
nionodisper.se  form  of  divided  metal  known  in  the  colloidal  si/e  regime.*  Sols  P  and 
Q  were  synthesised  by  alkaline  tetrakis  (hydromethyl)  phosphonittm  chloride  reduction 
of  chloroauric  acid/17)  All  the  hvdrosols  were  additionally  concentrated  by  removal  of 
the  solvent  (water)  before  measurement  of  the  nonlinear  optical  properties.  For  sols  P 
and  Q,  soduim  citrate  was  added  and  the  sols  boiled  for  15  minutes  before 
concentration.  For  all  the  sols.  PVP^  was  added,  to  a  concentration  of  about  2  a/mM 

of  metal  atoms,  before  boiling  down  using  a  hot  plate,  either  in  a  stream  of  nitrogen  or 
air. 


All  the  chemicals  used  were  analytical  reagent  grade,  except  tetrakis 
(hydromethyl)  phosphonium  chloride  which  was  only  a  general  purpose  reagent 
(Fluka).  Chloroauric  acid  was  used  as  obtained  from  Johnson  Matthey  and  contained 
49.419c  metal,  by  weight.  Poly(N-vinyl-2-pyrrolidonc)  at  40.000  molecular  weight 
(PVP40k)  and  at  10.0(K)  (PVP^J  were  used,  as  supplied  by  Polyscicnccs  inc..  All 

water  used  as  the  dispersion  medium  of  the  sols  was  distilled  twice,  latterly  from  an  all¬ 
glass  apparatus. 

2  2  MATERIAL  CHARACTERISATION 


The  materials  were  characterised  using  Transmission  Electron  Microscopy 
(TEM),  linear  UV-Visible  absorption  spectroscopy  and  a  range  of  nonlinear  optical 
techniques.  The  following  two  subsections  describe  the  linear  characterisation  methods. 
i.e.  TEM  and  (JV-Visiblc  absorption  spectroscopy. 

2.2.1  TRANSMISSION  ELECTRON  MICROSCOPY  : 

Transmission  electron  microscopy  is  an  extremely  useful  technique  for 
determining  the  size  and  shape  of  nanometric  particles  and  is  particularly  attractive  in 
this  case  because  of  the  high  contrast  afforded  by  the  interface  between  the  gold 
inclusions  and  the  surrounding  water.  Transmission  electron  microscopy  was 
performed  using  an  Hitachi  H7000  transmission  electron  microscope  operating  at  100 
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K  V  accciciaiion  voltage.  .Small  volumes  oi  ilia  c> 'iiccntiaicd  .so  is  wore  diluted  M  a  metal 
concentration  of  about  I  -  -  mM  with  doubly  distilled  water  and  a  drop  of  each  placed 
tor  several  seconds  in  contact  with  a  thin  non- perforated  carbon  film  supported  on  a 
v<)5  nun  copper  and.  The  particle  profiles  were  measured  according  to  their 
equivalent  circle  diameter",  and  the  resulting  distributions  of  diameters  assessed  for 
their  means  and  coefficients  of  variation  (standard  deviation  of  the  distribution 
expressed  as  a  percentage  fraction  of  die  mean).  ( mfortunaiely.  there  are  certain  errors 
a  uh  electron  microscopy  that  are  hard  to  estimate  : 

;  i  Inaccurate  magnification  by  the  microscope. 

in  Sampling  error. 

iii)  Artefact  of  size  introduced  bv  imaging  conditions  (extent  of  defocus,  si/e  of 
objective,  apertures  etc.)  and 

iv)  Measuring  inaccuracies  (particularly  induced  by  estimating  equivalent  circle 
diameters  of  aspherical  particles). 

5%  is  a  reasonable  value  for  the  error  in  the  absolute  value  of  <d>  (average  diameter) 
for  monodisperse  sols.  Since  most  of  the  errors  above  are  systematic  (apart  form 
sampling  errors  and  inhomogeneities  which  can  be  partially  estimated  by  examining 
different  micrographs  and  are  less  for  monodisperse  sols),  comparisons  between  one 
sol  and  another  still  hold.  For  polydisperse  sols  the  ahsoluie  error  is  greater. 

Figures  2.1  a-c  show  typical  micrographs  obtained.  As  is  evident,  the  si/e 
distribution  of  the  particles  is  very  narrow  and  they  can  be  considered  to  he  spherical  to 
a  high  degree  of  accuracy.  Micrographs  taken  of  the  same  samples  over  a  year  later 
showed  dial  they  were  free  of  aggregates  and  had  the  same  characteristics  as  the  freshly 
prepared  sols. 

2.2.2  UV-VISIBLE  ABSORPTION  SPECTROSCOPY  : 

1J  V- Visible  absorption  spectroscopy  is  a  standard  technique  for  obtaining  the 
energies  of  interband,  intraband  and  other  transitions  in  materials.  Literally  100's  of 
papers  have  been  published  on  the  linear  optical  properties  of  small  metal  particles.  See. 
for  example  references  3  and  4.  The  following  section  discusses,  therefore,  only  briefly 
the  general  features  of  die  linear  absorbance  spectrum  of  colloidal  gold. 

Hlcctrtc  dipole  absorption  in  Gold  composite  materials  results  in  a  peak  in  the 
visible  part  ol  the  absorption  spectrum,  called  the  "surface  plasmon  resonance"  peak 


y 


1 1..M  Micrographs 
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I  lk-  absorption  is  din.'  to  a  collective  electron  plasma  oscillation  (  piasmon  i  dial  ,s 
coupled  to  an  external  transverse  electromagnetic  Held  through  the  particle  surface.’ Is) 
The  piasmon  consists  of  the  conduction  electrons  of  the  metal  and  can  he  classically 
described  as  a  spherically  bounded  free  electron  wave  that  is  damped.  The  UV  belli 
absorbed  by  Silver  and  Gold  particles,  on  the  other  hand,  is  due  to  inierhand 
transitions,  and  is  quite  distinct  from  intrahand  excitation  of  the  piasmon  resonance  in 
the  visible.  The  fact  that  this  free  electron  part  of  the  optical  response  is  distinct  from 
the  bound  electron  and/or  interband  contribution  and  occurs  in  the  visible,  is  why  these 
colloids  exhibit  such  strong  colours. 

Depending  upon  the  size  of  the  inclusions,  two  different  theories  are  used  to 
explain  the  linear  absorption  spectra  of  the  composite  materials.  The  first  of  the  two.  the 
quasi-stalie  or  dipole  approximation  approach  is  valid  when  the  size  of  the  particle  is 
much  smaller  than  the  wavelength  of  light  incident  on  it.  This  reduces  the  problem  to 
that  of  a  panicles  immersed  in  a  uniform,  but  time  dependent,  field.  Thus  retardation 
can  be  neglected,  and  only  dipole  modes  need  to  be  considered.  The  quasi-static 
approximation  is  reasonable  for  particles  with  dimensions  <  Wnm  for  X  =  530nm. 

Within  this  approximation,  features  of  the  absorption  spectrum,  such  as 
position  and  width  of  the  surface  piasmon  peak  arc  size  independent.  There  is, 
however,  a  theory  that  the  dielectric  function  of  a  small  metal  sphere  is  dependent  upon 
its  size.  This  function  is  based  on  a  classical  model  of  a  reduced  electron  mean  free  path 
or  a  quantum-mechanical  model  that  leads  to  discrete  electron  energy  levels.  Both 
viewpoints  predict  that  the  imaginary  part  of  the  dielectric  function  will  increase  with 
decreasing  particle  size,  leading  to  a  broadening  of  the  surfacc-plasmon  width  for 
smaller  particles.  The  quantum  mechanical  model  also  predicts  a  slight  blue  shift  in  the 
peak  position  of  the  piasmon  resonance  with  decreasing  particle  size. 

When  the  size  of  the  inclusions  becomes  larger,  retardation  effects  have  to  be 
included.  At  this  stage  the  "Mic  theory"  <5>  is  used  to  explain  the  various  features  of  the 
absorption  spectrum.  According  to  the  Mie  theory,  the  piasmon  resonance  broadens 
ar.  j  the  peak  position  red  shifts  with  increasing  particle  sizes  -  similar  to  the  quantum 
mechanical  predictions  for  small  particles.  The  theory  derives  from  a  solution  of 
Maxwell’s  equations,  with  boundary  conditions  defined  such  that  the  electric  and 
magnetic  fields  just  outside  the  sphere  are  the  same  as  those  just  inside  ,  >.  the 
radiation  must  be  continuous  across  the  sphere  boundaries.  Using  these  boundary 
conditions,  it  is  possible  to  come  up  with  an  expression  relating  the  absorption 
coefficient,  a,  to  the  wavelength  of  the  exciting  radiation.  It  is  given  by. 


1 1 


u  = 


CD 


I.S.T  >' 

-7-1*,  - - - - - 

A  K  +  2M  +  2tv 


where  £,  and  f,  arc  the  real  and  imaginary  parus  of  the  dielectric  constant  of  the 
inclusions,  />  is  dieir  volume  f  raction  and  rh  is  the  (  purely  real  )  dielectric  constant  of 
die  host  dielectric. 


f  igures  2.2  and  2.3  show  the  linear  absorbance  spectrum  of  some  of  the  gold 
colloidal  samples  in  a  spectrophotometer  cell  of  path  length  I  mm. 


Wavelength  (nm) 


Figure  2.2  :  .inear  absorbance  spectrum  of  two  suniJarly  sized  gold  colloids 
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1'jgUTC  2.3  :  Linear  absorbance  spectrum  ol  two  gold  colloids  of  completely 
dilferent  si/es 


As  is  evident  from  the  figures,  it  is  difficult  to  come  to  any  definite  conclusion 
regarding  the  relationship  between  particle  size  and  width  of  the  absorption  band  and 
peak  position  of  the  plasmon  resonance.  It  is  probable  that  the  discrepancies  are  at  least 
partially  due  to  the  fact  that  most  theoretical  models  arc  based  on  idealised  samples  i.e. 
an  assembly  of  non-interacting  spheres  of  the  same  size.  In  reality,  of  course,  this  is 
not  the  case.  Unlike  silver  particles  which  have  a  sharp  surlacc-plasmon  resonance,  the 
resonance  in  gold  is  broad  because  of  the  relatively  large  bulk  value  of  the  imaginary 
part  of  the  dielectric  function,  e, In  addition,  there  is  the  onset  of  an  inlerband 

transition  in  gold  near  the  plasmon  resonance  that  leads  to  large  absorption  in  the  blue 
wing  of  the  absorbance  spectrum.  The  distribution  of  particle  sizes  and  shapes/^)  the 
number  of  lattice  defects (21)  and  the  large  value  of  the  bulk  dielectric  constant,  c ,  will 
tend  to  mask  changes  in  the  linear  absorbance  spectra  that  are  due  to  a  size  dependent 
£,  .The  host  matrix  also  plays  an  important  role  in  determining  the  width  and  the  peak 

position  of  the  surface  plasmon  peak. 

2.3  SIZE  EFFECTS 

Small  metal  particles  almost  never  have  plasmon  peaks  as  sharp  as  those 
predicted  by  the  simple  Mic  theory/22^  The  reason  for  this  is  that  the  fundamental 
electronic  properties  of  small  metal  particles  differ  from  those  of  the  bulk,  i.e.,  small 
metal  particles  arc  intrinsically  different  from  a  similar  volume  of  bulk  metal,  because  it 
is  the  near  infinite  continuous  lattice  of  metals  that  defines  their  electronic  structure. 
This  is  known  as  the  size  effect/4-  23-  24)  It  is  an  intrinsic  effect,  due  to  the  change  in 
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l Ik'  |Mi>|vHk\s  ol  ilk*  mnicnui  at  .small  crystallite  size.  whereas  i Ik*  \.h  i.ilum  in  optical 
pro|vi  ik's  with  paruele  diameter  (  as  described  by  the  Mie  tlieoi  v  )  is  an  i-.unn.sn •  effect 
defined  by  the  geometry'  of  the  interaction  of  light  with  a  sphere,  the  nature  of  the  metal 
remaining  unchanged. 

The  size  effect  al lecis  the  inclusions  in  two  distinct  ways,  firstly,  the 
delocalized  conduction  electrons  of  the  bulk  find  themselves,  in  the  inclusions, 
confined  to  regions  that  are  much  smaller  than  their  delocalization  lengths.  (Tor  an 
unconlined  metal,  delocalization  lengths  extend  over  several  unit  cells).  This  results  in 
discrete  optical  resonance  s  whose  position,  oscillator  strengths  and  dynamics  depend 
on  the  extension  of  the  artificial  confinement.  This  effect  is  known  as  the  Quantum 
confinement  effect  (  it  is  quantum  mechanical  in  nature  ).  So  basically,  it  relates  to  the 
alternations  in  the  energy  structure  of  the  material  and  comes  into  play  when  the 
physical  dimensions  of  the  material  become  comparable  with  the  characteristic  lengths 
that  govern  quantum  mechanical  processes. 

The  second  effect  that  stems  from  the  small  size  of  the  particles  is  that  the 
electric  field  that  acts  on  and  polarises  them  is  vastly  different  from  the  external 
macroscopic  electric  field  applied  to  the  composite  material.  This  is  due  not  only  to  the 
particle  size  but  also  because  their  dielectric  constant  is  completely  different  from  that  of 
the  surrounding  matrix.  This  effect  is  known  as  the  dielectric  confinement  effect 

Quantum  and  dielectric  confinement  effects  will  be  discussed  in  more  detail  in 
the  following  two  subsections.  They  can  be  used  to  enhance  the  nonlinear  response  of 
the  composite  material  and  are  particularly  sensitive  in  the  optical  frequency  range. 

2.3.1  QUANTUM  CONFINEMENT  : 

Quantum  confinement  in  small  metal  particles  was  discussed  for  the  first  time  in 
a  paper  by  Frdlich  published  in  1937/25)  He  showed  that  metallic  matter  in  the  form  of 
sufficiently  small  grains  behaves  qualitatively  differently  from  the  bulk  metal.  The  basis 
of  his  treatment  of  the  electronic  properties  of  these  particles  was  the  fact  that  the 
spacing  between  adjacent  energy  levels  increases  with  decreasing  particle  size.  He 
assumed  that  for  a  given  particle  size,  the  spacing  between  adjacent  levels  is  constant 
and  independent  of  energy.  The  theory  is  still  basically  unchanged  except  that  the 
artificial  concept  of  equally  spaced  energy  levels  is  removed  and  a  number  of  more 
realistic  assumptions  arc  introduced. 

Metal  nanocrystals  occupy  a  position  intermediate  between  a  molecule  and  the 
hulk  crystal.  Therefore,  the  choice  of  a  model  that  accounts  for  the  coexistence  of 
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le. tunes  Irom  (null  extremes  is  quite  complex.  I  or  a  system  consisting  ol  very  lew 
atoms,  the  most  usual  approach  is  the  one  using  molecular  orbitals,  hut  this  becomes 
enormously  complicated  as  the  number  of  atoms  is  increased.  On  the  other  hand,  the 
description  is  greatly  simplified  for  an  infinitely  extended  periodic  system  where  the 
behaviour  of  an  electron  can  be  determined  by  studying  the  appropriate  Sehrddinger 
equation  which  is  given  by 


HH '(/-)  = 


k  2  in 


VJ  +  V(r)  Y'(r)  =  EV{r) 


(2.2) 


where  V(r)  is  the  crystal  potential  experienced  by  the  electron  of  mass  in  .  and  */'(  r) 
and  E  are,  respectively,  the  state  function  and  energy  of  this  electron/^) 


Because  the  ions  in  a  perfect  crystal  are  arranged  in  a  regular  perfect  array,  the 
potential  V(r)  has  the  periodicity  of  the  underlying  Bravais  lattice,  i.e. 


V{r  +  R)=V{r)  (2.3) 

for  all  Bravais  lattice  vectors  R. 

The  problem  of  an  electron  in  a  solid  is,  in  principle,  a  many-clectron  problem 
for  the  full  Hamiltonian  of  the  solid  contains  not  only  the  one-electron  potentials 
describing  the  interactions  of  the  electron  with  the  atomic  nuclei,  but  also  pair  potentials 
describing  the  electron-electron  interactions.  In  the  independent  electron  approximation 
these  interactions  are  represented  by  an  effective  one-electron  potential  V(/). 

One  is  thus  led  to  examine  general  properties  of  the  Sehrddinger  equation  for  a 
single  electron,  i.e.,  equation  2.2. 

Independent  electrons,  each  of  which  obeys  a  one-electron  Sehrddinger 
equation  with  a  periodic  potential,  are  known  as  Bloch  electrons  W)  (  ^  opposed  to 
free  electrons  to  which  Bloch  electrons  reduce  when  the  periodic  potential  is  identically 
zero). 


According  to  the  "Bloch  theorem",  the  eigenstates  V7  of  the  one-electron 
Hamiltonian,  above,  can  be  chosen  to  have  the  form  of  a  plane  wave  times  a  function 
with  the  periodicity  of  the  Bravais  lattice  as  follows. 


^Jr)  =  c‘k'UJr) 


12. 4) 
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Wlii'ii'  t  '  .C  +  ft)  I  ',(/).  \vhcic  ft  is  any  l.iiiisv  wctor  I’lii*  coi  responding  ('and 
energy  l.r  ( k  )  is  a  reciprocal  space  periodic  (unction  L \{k  +  A  )  =  Fr(k )  where  K 

is  any  reciprocal  lattice  vector  and  k  is  the  wave  vector  that  labels  the  electron  state  in 
the  hand  /;  within  die  first  Brillouin  /one. 

l'or  an  infinite  crystal,  equation  2.4  has  the  form  of  a  tree  wave.  If,  however, 
the  periodicity  is  broken  te.g.  due  to  crystal  boundaries  or  defects  )  then  one  has  to 
represent  the  electron  motion  as  a  wave  packet  of  Bloch  states.  If  the  defect 
encompasses  several  unit  cells,  remains  essentially  unaffected,  and  close  to  k  =  0 

one  can  write. 


Y=r„(r)Un(r)  (25) 

Where  Fn(r)  is  the  envelope  function  that  has  to  be  solved  for. 

Metals  are  characterised  hv  a  single  half-filled  hand  up  to  the  Fermi  level  EF 
with  electron  and  hole  states  on  either  side  of  it.  The  electron  and  hole  wave  functions 
will  thus  he  of  the  form, 

=  ^t.r,,)(/f(re,rh)Ut(rt.rt)  (2>6) 

where  Ut  and  Uh  are  the  cell  periodic  part  of  the  Bloch  band  states  for  the  conduction 
(electron)  and  valence  (holes)  bands  respectively  at  k  =  0.  Equation  2.2  can  therefore 
be  rewritten  as. 


2m, 


r(rf.rh)  =  EF(r(,rh 


(2.7) 


Where,  V,h  =  Electron  -  hole  coulomb  interaction.  Vmh  =  residual  electron  and  hole 

coulomb  exchange  interactions,  which  arc  lumped  together  under  the  term  many  body 
effects,  and  Vt  =  Electron(hoJc)  -  impurity  interaction  potential. 


and  V,  collectively  constitute  V(r).  If.  in  addition  to  the  above 
potentials,  the  electrons  and  holes  find  themselves  confined  to  lengths  less  than  their 
delocalization  lengths,  (  c.g.,  in  crystallites  the  motion  of  the  charged  particles  is 
hindered  by  the  interface  with  the  dielectric  host  matrix  )  then  an  extra  potential 
describing  the  confinement  has  to  be  introduced  into  the  one-electron  Hamiltonian  as 
follows. 
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( 2.S ) 


-V.  +  V^  +  V^  f  V.  +  HM/'V..,/,.)  =/:T(r, 


Hie  confinement  potential  IT  can  he  visualised  as  an  infinite  spherical  potential  well. 
/.(■..  If  =  0  lor  r  <  a  and  If  -  ^  for  r  >  a  (  a  is  Lite  radius  ol  the  crystallite  ). 

flic  above  equation  can.  in  general,  only  he  solved  with  numerical  methods  and 
even  so  with  limited  usefulness  because  of  the  uncertainties  in  the  definitions  and 
determination  of  the  different  potential  terms.  However,  the  essential  features  of  the 
quantum  confinement  can  he  quantitatively  accounted  for  if  some  simple  analytical 
forms  for  these  potentials  are  introduced,  and  in  addition  their  relative  strengths  are 
properly  accounted  for. 

The  relative  impact  of  the  different  potential  terms  depends  on  the  extension  of 
the  confinement  with  respect  to  the  characteristic  lengths  that  are  associated  with  these 
potential  terms  and  determine  their  respective  strengths.  Whenever  the  confinement 
extension  is  larger  than  any  of  these  lengths,  the  energy  spectrum  is  determined  by  the 
corresponding  potential  term  while  the  confinement  only  acts  as  a  perturbation.  If  on 
the  other  hand,  the  confinement  extension  is  less  than  these  lengths,  the  roles  arc 
reversed  and  the  confinement  suppresses  the  different  potential  terms  and  imposes  its 
own  spectral  distribution/12’ 131 

In  the  case  of  metals,  it  turns  out  that  the  confinement  potential  W  completely 
dominates  all  other  potentials.  The  reason  for  this  is  that  the  half  filled  band  in  metals, 
which  is  usually  formed  with  s-  and  p-  orbitals,  can.  for  most  purposes,  be  replaced  by 
an  equivalent  pair  of  parabolic  bands,  situated  on  either  side  of  the  Fermi  level,  that  are 
mirror  images  of  each  other  and  touch  at  k  =  0.  The  wave  vector  dependent  dielectric 
constant  e(k)  is  infinite  at  k  =  0,  leading  to  the  potentials  Vrh,  Vmb ,  and  V,  being 
completely  screened  to  within  a  distance  rF  =  ,  which  is  of  the  order  of  a  lew 

angstroms,  i.e.  roughly  equal  to  the  lattice  constant.  Thus  electrons  and  holes  behave  as 
free  non-interacting  particles  over  any  distance  in  the  perfect  crystal. 

The  one-electron  periodic  Hamiltonian  can  therefore  be  replaced  by  the  free 
electron  and  hole  Hamiltonian,  leading  to  the  following  equation 


h2  t,2 

2  m..  2  m. 


F(rt.rh)  =  EF(rt,rt,) 


(2.9) 
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Ilk'  muon  I.'  this  iqu.idon  vields  a  signil icami v  ihllciviit  cigcnl nikUnii  and 
eigcuoncrgy  spectrum  from  ihc  hulk.  i.i’..  ihc  initially  continuous  energy  spectrum  ol 
die  tree  electron  Hamiltonian  is  replaced  by  a  discrete  energy  spectrum,  namely  one  of 
quantum  confined  levels  whose  spacing  depends  on  the  extent  of  the  confinement  and 
can  accordingly  he  modified  as  desired. 

In  spherical  co-ordinates,  equation  2.9  can  he  solved  to  give 


/•'(r.d.o) 


(2.10) 


=  RJ7{0.<t>)  (2.11) 

using  the  boundary  condition  F(r  =  a)  =  0. 

j,  is  the  spherical  bessel  function  of  order  /,  arJ  is  its  nth  zero  and  }',m(0,<j)  are  the 
spherical  harmonics.  Die  energy  of  the  state  given  by  quantum  numbers  n,  1  and  m  is. 


where  E0  = 


JL 

2  nra2 


is  the  ground  state  energy. 


2.3.2  DIELECTRIC  CONFINEMENT  : 


(2.12) 


Despite  their  disparity,  all  composite  materials  that  are  formed  by  uniformly 
dispersed  metal  or  semiconductor  crystallites  in  a  liquid  or  solid  transparent  dielectric 
share  an  important  feature  that  has  an  essential  impact  on  their  properties  in  the  optical 
frequency  range:  because  the  size  of  the  crystallites  is  much  smaller  than  the  incident 
wavelength  and  their  dielectric  constant  is  very  different  from  that  of  the  surrounding 
transparent  dielectric,  the  electric  field  that  acts  on  and  polarises  the  charges  of  these 
crystallites  can  be  completely  different  from  the  macroscopic  Maxwell  field.  This  effect 
is  known  as  the  dielectric  confinement  effect 

2.3.2. 1  LINEAR  REGIME  : 

In  general,  the  metallic  crystallites  embedded  in  the  dielectric  arc  not  completely 
uniform  either  in  size  or  in  shape.  However,  in  order  to  explain  the  essential  features  of 
dielectric  confinement  several  features  of  the  composite  material  can  be  considerably 
simplified,  firstly,  it  can  he  assumed  that  the  surrounding  host  dielectric,  liquid  or 
solid,  is  an  ideal  isotropic  dielectric  of  dielectric  constant  ch  with  the  metal  particles 
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umli'inilv  and  randomly  distributed  in  .small  volume  concentration  in  u.  They  are 
assumed  to  he  spherical  in  shape,  with  a  diameter  d  =  2a  that  is  much  smaller  than  the 
optical  wavelength  X. 


Figure  2.4  :  A  small  sphere  of  dielectric  constant  f  embedded  in  a  matrix  of 
dielectric  constant  f, . 

In  the  linear  regime  the  relevant  optical  coefficient  of  such  a  crystallite  or 
inclusion,  of  volume  V,  is  the  polarizability  av  and  its  dielectric  constant  £(  can  he 

defined  by  the  relation. 

£t .  =  1  +  47T£(<I)  =  1  +  4n~  (2.13) 

where  is  the  linear  susceptibility. 

The  dielectric  constant  is  expected  to  be  a  function  of  the  particle  size.  Its  limit  for  large 
panicles  is  the  bulk  value  cR . 

Secondly,  if  it  is  assumed  that  the  volume  density  of  the  particles  in  the 
dielectric,  p ,  is  so  small  that  the  interaction  of  the  induced  dipole  moments  on  the 
particles  can  be  ignored,  then  an  effective  dielectric  constant  c  for  the  composite 
medium  can  be  introduced,  whose  relation  to  eh,  and  p  is  given  by  the  Maxwell- 

Garnet  expression/6) 


c-e»  £,-c> 

£-  +  2£„  £,  +  2  eh 


(2.14) 


This  relation  is  a  straightforward  consequence  of  the  Clausius-Mossotti 
approximation  for  the  local  field  corrections  for  spherical  polarizable  particles  (26f  For 
/>  «  1  it  reduces  to. 
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r  -  r,  +■  3/uy 


v C  +  ~Ch 


(2.15) 


liquation  2.14  can  he  derived  from  Lorent/.-Lorenz  local  lield  arguments  as 
follows: 

Consider  a  collection  of  molecules  subjected  to  an  external  electric  field,  E.  The 
dipole  moment  induced  in  any  one  ol' the  molecules  is  given  bv. 


/'  = 


(2.16) 


where  Elocal  is  the  Lorentz  local  field,  which  consists  of  the  external  field  and  the  Held 
due  to  the  other  molecules,  and  a  is  the  linear  polarizability. 


4,*  =  E+yP 


and  P  is  the  induced  polarisation. 


(2.17) 


If  the  number  of  molecules  =  N ,  then. 


—  4  7T 

P-Nji-  NaElocal  =  Na\  E  +  —  P 


(2.18) 


or 


P  = 


Na 


1 - Na 

3 


E  =  xwE 


(2.19) 


where. 


*U)  = 


Na 


ATC  X, 

- Na 

3 


Equation  2.20  can  be  rewritten  as, 


£(l)-l  \K 

—r - =  — Na 

£n)  +  2  3 


where  c  '  is  defined  as. 


(2.20) 


(2.21) 
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£"'  -1+4 K'/_ 

filiation  2.21  is  the  Lorentz-Lorenz  law.  It  can  alternatively  he  stated  as. 


.  rl  '  t-  2 

/  '  =  -  Na 

3 


(2.24) 


e"1  +2 


(2.24) 


f  £(l)  +  2  | 

where  -  is  called  die  local  field  enhancement  tactor. 

3 


These  arguments  can  also  he  applied  to  a  composite  material  in  the  presence  of  a 
uniform  electric  field  29>  The  polarisation  induced  in  the  material  due  to  the 
external  E  field  is  given  by. 


P^XhE+P 


(2.25) 


where  Xh  is  susceptibility  of  the  host  matrix  and  P  is  the  additional  polarisation 
due  to  the  inclusions. 


P  =  Nfi  =  NehaE, 


(2.26) 


^  is  the  dipole  moment  of  the  inclusions,  a,  =  a'  — — —  ,  is  their  polarizability  and 

r  4- 1  r 


£.  +  l£. 


E ,  is  the  local  field,  i.e..  the  field  in  their  vicinity. 


Just  as  in  the  derivation  of  the  Lorentz-Lorcnz  law,  one  can  write. 


E,  =  E  +  —  P 
3ck 


(2.27) 


which  implies  that. 


I - —  /Va, 
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(2.28) 


Combining  equations  2.25  and  2.28  one  gets. 
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I isiiiii  the  definition.  <•  =  i  +  4^,  tor  the  dielectric  constant,  equation  2.30  can  he 
alternatively  expressed  as  : 


£  +  2t',. 


\  £,  +  2c 


ft  / 


(2.31) 


which  is  the  same  as  equation  2. 14. 

£,  is  complex.  c,  =  c ,  +  /£, ,  and  frequency  dependent.  It  can  he  seen  that  the 
expression  in  equation  2.31  shows  an  enhancement  close  to  the  frequency  at. ,  such 
that. 


£,((0,)  +  2£*=0  (2.32) 

which  is  the  condition  for  the  surface  excitation  or  surface  plasmon  resonance 
frequency.  The  width  of  this  resonance  is  determined  by  £,  .The  absorption  coefficient 

as  a  function  of  wavelength  can  easily  he  calculated  from  equation  2.3 1 .  It  is  given  by. 
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(2.33) 


which  is  the  same  as  the  result  that  is  obtained  using  the  Mie  theory  to  calculate  the 
optical  constants  of  composite  materials. 


Combining  equations  2.27,  2.28  and  2.31  it  is  possible  to  derive  an  explicit 
relation  between  the  local  field  experienced  by  the  inclusions  and  the  external 
macroscopic  which  is  given  by. 


£,= 


3£, 


Similarly,  the  field  inside  the  inclusions  is  related  to  the  external  field  through. 


-n 


(2.34) 


(2.35) 


Combining  equations  2.34  and  2.35  the  Held  inside  the  inclusions  can  he  expressed  in 
terms  ol  the  local  field  as. 


^l  =  £LT=£L^=>E„=flE, 


(2.36) 


where  /,  is  called  the  local  field  enhancement  factor,  liquation  2.36  constitutes  the 
main  result  of  this  section. 


2. 3. 2.2  NONLINEAR  REGIME  : 


In  the  presence  of  an  intense  electric  field  the  polarisation  induced  in  the 
crystallites  may  be  written  as. 


AP  =  P0)  +  P{2)  +  P0) (2.37) 

where  Pt'n)  with  n>  1  is  the  nonlinear  polarisation  term  of  order  n.  In  the  case  of  an 
isotropic  composite  with  random  distribution  of  inclusions  P{2n)  =  0, 


P(l)  =  x0)E,  and 
p(3)  =XWEEE, 

where  2(1)  and  2f<3)  are  the  linear  and  third  order  effective  susceptibilities  of  the 
composite  material. 

The  term  Pi^  is  related  to  third  order  effects,  the  most  important  of  which  is  the 
intensity-dependent  refractive  index  or  the  optical  Kerr  effect.  This  effect  is  involved  in 
a  wide  variety  of  important  processes,  including  self-focusing  of  a  laser  beam/313  self¬ 
phase  and  frequency  modulation/32)  "soliton"  pulse  propagation/33)  and  "optical  phase 
conjugation"/343  It  can  be  described  as  an  optically  induced  change  of  the  optical 
dielectric  constant,  i.e. 


8e  =  \2kx0)\E{o))\  (2.38) 

This  change  of  £  contains  contributions  from  both  the  surrounding  dielectric  and  the 
embedded  crystallites,  denoted  by  8t:h  and  Se,  respectively.  Close  to  the  surface 

plasmon  resonance  frequency  the  contribution  from  the  surrounding  dielectric  is 
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negligible  compared  to  that  I'rom  the  crvv.,ihie.  even  lor  />•  ■  I.  Tims  a  change.  Si '  ol 
the  dielectric  constant  ol  ihc  inclusions  will  result  in  a  change  in  r.  the  dielectric 
constant  of  the  composite  given  by. 


Sc  =  i> 


V, 


e,  +  2c 


Sc,  =  p.t\'Sc 


h ) 


The  above  equation  is  obtained  by  differentiating  equation  2. 1 5. 


(2.39) 


If  the  third  order  susceptibility  relevant  to  the  internal  field  within  inclusions  is 
designated  by  _£(<3)lhen  by  analogy  to  equation  2.3N. 


Substituting  for  Sc,  in  equation  2.39  gives 


(2.40) 


8e  =  pfi15e,  =  [2pn\f]\-fcx':)n- 


(2.41) 


and. 


*(3,  =  pf/,|V,V  (2-42) 

Thus  it  can  be  seen  dial  the  nonlinearity  can  be  enhanced  by  an  amount  equal  to 
the  fourth  power  of  the  local  field  factor,  close  to  the  surface  plastnon  resonance 
frequency.  As  one  moves  away  from  this  resonance  frequency,  the  local  field  factor 
approaches  unity.  At  this  stage  dielectric  confinement  plays  no  role  in  the  enhancement 
of  the  nonlinearity. 

Another  important  point  to  be  noted  in  equation  2.42  is  that  the  macroscopic 
susceptibility  x0)  can  also  be  enhanced  by  the  intrinsic  quantum  confinement  mediated 
enhancement  of  the  susceptibility  of  the  inclusions,  x}3). 
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2.4  NONUNEAR  OPTICAL  PROPERTIES 


Hie  liiird  order  susceptibility  of  die  inclusions,  xl''  > oan  he  calculated  using  the 
quantum  mechanical  expression  in  equauon  1.2  which  is  repeated  here  for  convenience. 

t.Q.J  .‘I  J 

(  7lu  7Za 7I11 7tr  \ 

l  - - - - -  +  47  similar  terms 

{{^  -  <a,  -  -  w3)(A,,  -  O),  -  I  (2.43) 


In  order  to  calculate  the  susceptibility,  the  matrix  elements  of  the  transition  dipole 
moment,  n,  must  be  known.  They  can  be  calculated  using  the  relation 

[W,r]  =  -/j"  — where  p  is  the  momentum  operator  and  n  =  jp 

The  radial  part  of  r  between  two  states  r  (quantum  numbers  n,  1)  and  s  (quantum 
numbers  n',  1+1)  is  given  by, 

a 

f  r2drRnlrRH.M  (2.44) 

o 

Substituting  for  R ^  from  equation  2.11,  changing  the  integration  variable  from  r  to 

—  =  x  and  integrating  by  parts  gives, 
a 


(a2,  ~«L,) 


a  =  4  aE0 


{^nl*\  ~  E«l) 


Therefore  the  radial  pan  of  the  n  matrix  element  Is 


(2.45) 


4ae£0  {ErE,)V2 
ih  Es  -  Er 


(2.46) 


With  these  transition  dipole  moment  elements  it  is  possible  to  calculate  the  third  order 
susceptibility,  £<3\  of  the  composite  materials  using  the  expression  in  equation  2.43.  If 
the  experiments  are  performed  in  the  vicinity  of  a  resonance,  only  a  few  of  the  elements 
are  significant 
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in.si  contribution.  called  the  mirahand  contribution.  i.s  caused  hv  transitions  within  the 
s-p  conduction  hand.  It  is  mainly  imaginary  and  negative  and  shows  a  s trout:  si/e 
dependence.  I'll  is  is  to  he  expected  since  the  highly  delocalised  electrons  in  the 
conduction  hand  would  Iv  greatly  affected  hv  a  decrease  in  their  delocalisation  lengths 
due  tit  quantum  confinement.  The  second  contribution,  called  the  interhand 
contrihution.  is  due  to  transitions  Irom  the  d  -  valence  hand  to  the  s-p  conduction  hand. 
It  is  also  mainly  imaginary  and  negative  hut  is  independent  of  the  si/e  of  the  inclusions. 
There  are  two  main  reasons  lor  this,  f  irstly,  the  d  -  hand  electrons  have  a  much  greater 
effective  mass  than  free  electrons.  Secondly  they  tire  already  localised  and  are 
consequently  unaffected  hy  the  quantum  conlinement.  The  last  contribution,  called  the 
hot  electron  contribution,  arises  due  to  a  change  in  the  Fermi-Dirac  distribution  of 
electrons,  about  the  Fermi  level,  caused  by  their  increase  in  temperature  on  absorbing 
energy  in  a  resonant  process.  This  contribution  is  imaginary,  like  the  first  two,  but  it  is 
positive.  In  the  chapters  following  this  one.  several  experiments  are  described  which 
were  undertaken  in  order  evaluate  the  magnitude  ol  the  different  contributions  to  z'3> 

and  to  ascertain  which  one.  if  any,  was  the  dominant  one. 

The  next  three  sections  of  this  chapter  are  devoted  to  a  detailed  discussion  of  the 
three  different  contributions. 

2.4.1  THE  INTRABAND  CONTRIBUTION  : 

Equation  2.43  for  the  calculation  of  in  terms  of  the  transition  dipole 
moment  matrix  elements  contains  48  terms.  Keeping  only  those  terms  for  which  one  of 
the  factors  in  the  denominator  is  equal  to  /  /  Tj .  where  T,  is  the  lifetime  of  the  states, 
the  number  of  terms  reduces  to  32;  16  different  terms,  each  appearing  twice. 

Considering  the  different  components  of  the  ^lJ)  tensor,  the  terms  can  be 
divided  into  three  groups  :  eight  terms  in  four  terms  in  ;^3),  and  four  terms  in 
such  that. 


Til  =2(Tl3'+Tr)  +  Tl3’) 

y  1  -')  —  ")/  r*  y^  ^  4.  p  y  ^  ^  _l_  *%  y'  \ 
A  tvn  \  rvrvA  I  r  uwA  2  *  mi  A  3  J 


(2.47a) 


!.47b) 


Til  =  Til  =  *'™(t1  +  TVm)  +  ^  Lrvv  ( Ti 3 '  +  T\”)  +  ^,v„(t;3'  +  T.l3’)  (2.47c) 

where  i.s  a  coirective  multiplicative  factor  when  Ujkl)  t-  (  vvxv)  and 
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- Tr - i - T  f  v  i  ' ; 

+  (0  +  /TV1 )  +  ru)‘  +  7T' 

where  V  is  the  volume  of  the  crystallite,  7,  is  the  dephasing  lime,  A  is  an  angular 
form  factor,  n„  is  defined  by  equation  2.46  and  co!r  =  (£f  -  Er),  setting  ti  =  1  for 

convenience. 

There  are  two  main  contributions  to  The  first  one  arises  when 
O)  ~  Es  —  Er,  and  the  second  one  when  Es  ~  Er.  In  the  first  case  it  is  only  necessary  to 
consider  the  first  term  in  equation  2.49  since  its  contribution  dominates  that  of  the 
second  term.  In  the  second  case,  both  terms  need  to  be  considered.  However,  it  turns 
out  that  the  overall  contribution  to  £,<3)  in  this  case  is  much  less  than  the  in  first  case 
and  it  can  be  safely  ignored. 

Substituting  for  Kn  from  equation  2.46,  equation  2.49  can  be  rewritten  as. 


7*(Wl 

•  i  CO  r.s 


(E,-E,y 


for  the  case  when  Q)sr  =  (£,  -  Er).  The  summation  in  the  above  equation  can  be 
replaced  by  an  integral  to  give. 


(2.50) 
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where  v(li)  is  ihe  densiiv  of  suites.  In  the  second  integral.  the  densitv  - takes  into 

account  die  fact  dial  the  /  iiuinher  ol  the  s  state  is  determined  hy  the  /  tuitnher  ol  the  r 
slate  and  is  ol  different  parity. 


Equation  2.5 1  can  be  solved  to  give. 


X\»  = 

9  n  (0 


(2.52) 


with  v  =  —  and  jf.(v)  =  —  f ,r5/‘(.r  +  v)v7dx . 

r  1  it  * 


The  quantum  mechanical  expressions  for  and  X?]  arc  given  by. 


/,31  =  l\/Zi/f  X-  ~ 

-  5  t*  A 


3  T:  o)*  -  (o  -  f'771  )j(w(,  -  (of  +  T:2 
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((our  +  (o  +  iT;')\((o,r-(o)2  +  t;7 


(2.54) 


The  only  important  contribution  to  ^3)  is  when  (0ls  and  w  ,  are  close  lo  (0.  In  this 

case  u  is  only  necessary  to  consider  the  first  term  of  equation  2.53,  and  hy  proceeding 
in  the  same  way  as  lor  x\:'\  the  final  result  obtained  is. 


Xi''  =  i 
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(2.55) 


with 


Xi(v)-  J  r  '(*  +  v)  ''(x  +  2v)'n (Lx  (2.56) 

I- V 

Similarly,  lor  x\\  the  only  important  contribution  is  when  a)sr  and  (o n  are  close  to 
(O,  which  means  that  the  first  term  in  equation  2.54  can  be  ignored,  giving. 


X™  =  i 
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■AT, 


(4ae)4 
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-3/2 
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8i(v) 


(2.57) 


with 


*.(v)*-t}*M0,  +  v),n(*-v)"d» 
where  ar0  =  max(v,  I  -  V'). 


(2.58) 


Combining  equations  2.47a,  2.52,  2.55  and  2.57  yields. 


xZ  =  2  Xl3) 


I--2- 


(  ...  3 


1/2 


T 


2  £, 


2  V  /  ) 


(2.59) 


(2.60) 


where. 


a 


o 


g.(v) 

?2(V')  +  ^(V') 


or.  substituting  lor  jj1’. 


(2.61) 
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ill  SI  units. 

When  n  -  -*<». 

When  rt«n().  z'j’t  «<i 
When  a»a0,  **  a 

For  the  ease  of  gold,  with  Ef  =  5.5ev  and  setting  T2  =  20fs,  a0  ~  136  A  . 

The  magnitude  of  y' l|  can  he  approximated  by  inserting  suitable  values  for  the 
physical  parameters  in  equation  2.63  above.  Setting  7,  =0.5ps  and  (0  =  3.55x10’ V 
yields,  |*(3)|  =  4.23xlO',Vl'"2. 

2.4.2  THE  INTERBAND  CONTRIBUTION  : 

2.4.2. 1  LINEAR  REGIME  : 


It  is  well  known  that  the  generation  of  an  optical  plasma,  due  to  interband 
transitions,  can  result  in  a  strong  modification  of  the  dielectric  function  e(co)  in  a 
semiconductor.*36)  The  same  argument  can  be  applied  to  the  metallic  crystallites  in  the 
composite  material.  According  to  the  Drude  theory  of  metals.'37-  3S)  the  complex 
dielectric  constant  of  a  free  electron  metal,  to  a  first  approximation,  is  given  by 


(On 

e{co)  =  1 - r 

(O' 


(2.64) 


where  cop  = 


2  Ame1  . 


m 


is  the  plasma  frequency,  with  n,  e  and  m  being  the  electronic 


concentration,  charge  and  mass  respectively. 


In  the  case  of  alkali  metals,  there  is  a  striking  correspondence  between  theory 
and  experiment,  in  this  respect.  In  other  metals,  however,  different  contributions  to  the 
dielectric  constant  compete  quite  substantially  with  the  Drude  term. 
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The  dielectric  constant  of  noble  metals  consists  ol  two  terms:  one  due  to 
intraband  transitions  within  the  s  -  p  conduction  hand  and  the  other  due  to  interhand 
transitions  from  die  d  -  valence  hand  to  die  s  -  p  conduction  hand. 

(ot)  =  (2.05) 

The  free  electron  intraband  term  corresponds  to  die  Drude  dieleetnc  constant,  />.. 


1  Mr.,  '  f  ..  \ 

(oyco  +  il  r.„ ) 

where  Teff  is  the  mean  collision  time  for  free  electrons  taking  the  surface  into  account. 
In  the  large  sphere  limit  zrff  reduces  to  T2  /  2. 

The  imaginary  part  of  £,nler  is  dependent  on  the  joint  density  of  states  involved 
m  the  transition  and  on  the  momentum  operator  between  the  states.  It  is  given  by, 

(2.67) 

Ptj(i t) is  the  matrix  clement  of  the  momentum  operator  between  states  corresponding  the 
quasi  momentum  hk,  /,  and  fi  are  the  occupation  numbers  and  E{  and  £;  arc  the 

energies  of  the  initial  and  final  levels.  Near  the  surface  plasma  resonance  frequency,  die 
sum  over  the  indices  i  and  j  is  dominated  by  one  term  whose  initial  state  is  the  top  of 
the  valence  band  and  final  state  is  the  conduction  band.  Since  the  integral  in  k  space  is 
dominated  by  the  vicinity  of  the  X  point  of  the  Brillouin  zone  where  the  band  gap  is 
1.7eV,  (See  figure  2.5)  the  dependence  of  T(A.  )  is  weak  and  it  can  be  approximated  to 

be  a  constant  P.  Equation  2.67  can  then  be  rewritten  as. 


(2.68) 


and  the  corresponding  absorpdon  coefficient  a{(o)  is  given  by. 


,  ,  0) 

«(<»)* -7 TT-^nur  (2.69) 

n(0))c 

J(o))  is  the  joint  density  of  states.  In  thermal  equilibrium  the  contribution  to  the 
electronic  density  from  electrons  in  the  nth  band  with  wave  vectors  in  the  infinitesimal 
volume  clement  dk  of  k  space  is  given  by  the  Permi  distribution 
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Band  structure  of  cold  (lie.  2.5)  Reference  2.39 
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f(T)  = 


I 


(2.70) 


i  ,  *cr 
\+e  * 

with  SI:  =  l:h  +  hu)),  where  /:,  is  the  d-band  energy.  kK  is  Boltzmanns  constant 
and  T  is  die  temperature. 


Hie  joint  density  of  states  is  therefore  given  hy. 


J(0>)  =  1  j  (l  -  J  ,(k))S(l::  (k  )  -  I- (k )  -  hto) 

I  I  71 


(2.71) 


NONLINEAR  REGIME  : 

If  the  excitation  frequency  used  is  resonant  with  a  particular  interband 
transition,  the  main  contribution  to  %0)  will  consist  of  the  contribution  of  the  two  levels 
involved  in  that  transition.  Since  it  is  resonant,  it  will  be  mostly  imaginary  with  the 
imaginary  part  given  by  ??? 


(2.72) 


Where  T[  and  T2  are  the  energy  lifetime  and  dephasing  lime  corresponding  to  the  two 
level  system  and  A  is  an  angular  form  factor  =1/5. 

Since  only  contributions  from  the  X  point  of  the  Brillouin  zone  are  relevant  to 
the  integral  in  equation  2.72,  by  analogy  to  equation  2.68,  equation  2.72  can  be 
rewritten  as. 


■mja—j-^Tip^rWyC, 


(2.77) 


The  product  J(o)^P\2  can  be  obtained  from  the  value  of  e'inur  which  in  turn  can  be 
obtained  by  subtracting  the  Drude  term  from  the  overall  dielectric  constant  l^l2  is  given 
by  the  expression  W 


7  m 


(2.74) 
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wiK'iv  1 1  is  nil'  modulus  < *i  i hi*  reciprocal  iuiiicc  w ,t v i-  wctoi  iclcvaiii  me  ii.msiuon 
and  /a.-;  i.s  the  gap  energy.  faking  i'uM,  -  Id.  <»’  ~  I. '",!<>  \  //•  ( I  rum  i  Ho  point  X  in 

the  Briliouin  /one)  and  co  =  A55xl()‘ \v  1  (surface  plasmon  resonance  frequency), 
gives 


\n{  =  \l<)4,(Ag  in  v  )' 


(2.75 


and 


lm  -  — l.SxlO  'in  1'  ‘  t2.7o) 

In  order  to  obtain  the  result  in  equation  2.76.  the  values  lor  /',  and  7',  were 
approximated  to  be  2xl()'n  and  2x K)'1 4  respectively/^ 

This  contribution  to  x0)  *s  size  independent  down  to  very  small  sizes.  The 
reason  lor  this  is  that  the  quantum  si/e  effect  should  be  unobservable  for  the  d  -  levels 
because  of  their  large  effective  masses  and  the  fact  that  the  electrons  already  have  a 
localised  character.  The  calculation  of  the  interband  contribution  is  affected  by  s-p  band 
only  for  very  small  spheres  with  sizes  less  than  2.5nm/41-  42> 

Thus,  in  this  case,  the  contribution  to  xW  is  approximately  an  order  of 
magnitude  greater  than  that  due  to  intraband  transitions. 

2.4.3  THE  HOT  ELECTRON  CONTRIBUTION  : 

The  third  mechanism  contributing  to  the  susceptibility  of  the  inclusions  is  the 
hot  electron  contribution.  It  results  from  the  modification  of  the  population  of  the 
electron  states,  the  Fermi-Dirac  distribution,  caused  by  the  elevation  of  their 
temperatures  subsequent  to  the  absorption  of  photons  in  the  resonant  process,  but 
before  the  heat  is  released  to  the  lattice  of  the  crystallite;  this  leads  to  a  contribution  to 
X'})  that  is  positive,  imaginary  and  size  independent/9- ,0- I2- 13) 

Because  the  specific  heat  of  the  conduction  electrons  is  weak,  on  the  absorption 
of  photons  they  can  heat  up  to  a  temperature  that  is  much  higher  than  that  of  the  lattice. 
It  has  been  shown  that  it  takes  a  lew  ps  for  these  hot  electrons  to  come  back  into 
thermal  equilibrium  with  the  lattice/4-*-  44 )  During  this  lime,  their  Fermi-Dirac 
distribution  gets  modified  resulting  in  part  of  the  one-electron  levels  below  the  Fermi 
level  being  emptied  and  part  of  the  one-electron  levels  above  the  Fermi  level  being 
occupied.  This  leads  to  a  modification  of  the  dielectric  constant  of  the  inclusions  and 
consequently  contributes  to  the  overall  susceptibility. 


/4 


I  hi.'  cnerg\  iii  iIk'  photons  in  the  vicimtv  ol  the  surlace  plasmon  resonance 

frequency  is  very  close  10  the  hand  gap  which  exists  in  gold  near  the  L  point  ol  the 
Brillouin  /one  1 I  ho  .  =  2.AeV)  and  therefore  this  effect  is  very  strong.  The  relevant  pan 

of  the  total  dielectric  function  is  called  r,  and  its  imaginary  pari  is  given  bv,(1(h 

r]  =  C(\ -  f(T))  (2.77) 

where  f(T)  is  the  occupation  number  of  the  conduction  electrons  and  C  is  a 
temperature  independent  constant. 

The  change  in  the  dielectric  function  as  a  function  of  temperature  is  obtained  by 
differentiating  equation  2.77  : 


(2.82) 


uheic  ;  is  the  local  held  factor.  is  me  mi.igm.iiv  p.m  ol  ilu*  Unulc  contribution  i<> 
ihc  dielectric  constant.  r„  is  the  lime  constant  lor  the  cooling  of  the  electrons,  /,,  is  the 
incident  laser  energy  and  n  is  the  refractive  index  ol  the  composite  material. 

flic  quantity  rn  r„  inequation  2.K2  is  given  In. 


r„ 

(»'  r,, 


(2. Set 


where  ( o ,  is  the  plasma  frequency  and  the  ratio  is  equal  to  the  numher  ol  collisions 
'  '  r  „ 

an  electron  has  to  undergo  in  order  to  transler  its  excess  energy,  either  to  the  metallic 
particle  phonons  or  the  phonons  of  the  surrounding  dielectric  through  the  interlace.  If 
the  excess  energy  is  2.33  eV  ( luo)  and  each  collision,  on  average,  transfers  0.02eV. 
then  the  total  numher  of  collisions  required  is  approximately  12().<43)  Taking 
(0p  -  3. 49x10' V1  gives  r,,Tn  -  f3.\l()  "V 

Combining  equations  2.79,  2. SO.  2. SI  and  2.S2  gives  an  expression  for  the 
change  in  the  dielectric  constant  of  die  inclusions  which  is  related  to  their  susceptibility: 


Se, 


(0  24;r 
ncyr  2 4n: 


24  ir 
nc 


(2.84) 

(2.85) 


where 


yO)  _ 

x  dT  yr  24  k: 

Using  y  -  66 gives 


(2.8ft) 


X0>  -  10'V.vu  =  1.5xl0'lom'  V/  ;  (2.87) 

depending  on  the  wavelength  of  the  exciting  radiation. 

The  contribution  is  mainly  imaginary.  Contrary  to  the  first  two,  its  sign  is  positive. 
Also,  it  is  an  order  of  magnitude  greater  than  the  interhand  contribution  which  is.  in 
(urn.  an  order  of  magnitude  greater  than  the  intraband  contribution. 


CONCLUSIONS 

From  this  chapter  the  obvious  conclusion  would  be  that  in  an  experimental 
situation,  the  hot  electron  contribution  would  be  the  dominant  one  followed  by  the 
inierhand  and  iniraband  contributions,  in  that  order. 

This,  however,  is  not  strictly  true.  The  numbers  obtained  are  only  approximate 
as  a  lot  of  the  parameters  in  the  equations  are  not  known  to  a  sufficient  degree  of 
accuracy  and  do  not  merit  such  a  statement.  It  should  be  possible,  though,  to  determine 
experimentally,  from  the  sign  and  the  response  time  of  the  susceptibility,  the  physical 
mechanism  responsible  for  the  process. 

The  following  four  chapters  describe  different  experimental  techniques  for 
determining  £(3)  and,  hopefully,  the  physics  underlying  it. 
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CHAPTER  < 

intensity  oi:pi:ndi:nt  absorption  measurements 


Various  experimental  techniques  are  known  lor  determining  the  real  and 
imaginary  parts  of  the  third  order  susceptibility  x''! •  discussed  in  chapter  1.  This 
chapter  focuses  on  a  method  lor  evaluating  these  quantities  in  a  nonlinear  material  by 
measuring,  as  a  function  of  the  incident  light,  the  material's  absorption  and  the  phase 
change  induced  by  it/1-  -■  A  completely  linear  material  will  show  no  change  in 
absorption  with  intensity;  materials  that  are  "two  photon  (  or  multiphoton  )  absorbers" 
show  an  increase  in  absorption  with  increasing  intensity  and  materials  that  arc 
"saturable  absorbers"  show  a  decrease  in  absorption  with  increasing  intensity. 

EXPERIMENTAL  TECHNIQUE 

The  experimental  set-up  is  shown  in  figure  3. 1 .  Using  a  Gaussian  laser  beam  in 
light-focus  geometry,  the  transmittance  of  the  nonlinear  medium  is  measured,  in  the  far 
field,  as  a  function  of  its  position  Z,  measured  with  respect  to  the  focal  plane  of  the 
lens.  Such  a  trace  is  expected  to  be  symmetric  with  respect  to  the  focus,  where  there  is 
either  a  minimum  transmittance  (e.g.  multiphoton  absorption)  or  a  maximum 
transmittance  (e.g.  saturation  of  absorption).  The  coefficients  of  nonlinear  absorption 
and  hence  the  imaginary  pan  of  %0]  can  be  easily  calculated  from  such  a  curve. 

If  an  aperture  is  placed  in  front  of  the  photodiode,  it  is  possible  to  detect 
contributions  from  the  real  part  of  the  susceptibility  as  well.  This  contribution  causes 
the  beam  to  he  focused  or  defocused  (  depending  on  the  sign  of  the  nonlinearity  )  as  it 
passes  through  the  sample.  Focusing  of  the  beam,  when  the  sample  is  on  the  positive 
side  of  the  focal  plane,  causes  increased  aperture  transmittance,  while  defocusing  has 
the  converse  effect.  If  the  sample  is  on  the  negative  side  of  the  focal  plane,  the  situation 
is  reversed.  By  subtracting  the  effect  of  the  imaginary  component  from  this  curve  (  by 
using  the  first  trace  )  it  is  possible  to  separately  evaluate  both  the  components  of  the 
nonlinear  susceptibility/4- 5* 

The  laser  system  used  was  a  tuneable  PRA  Nitrogen  pumped  dye  laser  with 
5(X)ps  pulses  operating  at  a  repetition  rate  ol  5  Hz.  The  experiments  were  performed  at 
two  different  wavelengths  of  516nm  and  522nm,  which  lie  very  close  to  the  surface 
plasmon  resonance  frequency.  The  dye  used  was  Coumarin  4X5.  The  output  energy 
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from  ihc  laser  at  these  wavelengths  was  -  40//./  and  the  beam  was  focused  down  to  a 
spot  size  of  -  40/iifj  giving  a  peak  intensity  of  the  order  of  lOMf/ cm*. 

Focussing  Lens 


Motorised  Translation  Stage 


Figure  3.1  :  Experimental  Set-up  for  Intensity  Dependent  Absorption 
Measurements 


GENERAL  THEORY 

The  description  of  resonant  optical  processes  is  greatly  simplified  by  restricting 
attention  to  the  dominant  resonant  transition.  The  most  widely  used  description  is  based 
on  a  model  of  a  two  level  atom.  If  there  is  a  distribution  of  transition  frequencies 
between  the  levels  the  system  is  considered  to  be  inhomogenously  broadened. 
Saturated  absorption  spectroscopy  is  a  useful  tool  for  determining  whether  a  material 
can  be  described  by  homogenous  or  inhomogenous  saturation  because  the  dependence 
of  the  absorption  coefficient  on  intensity  is  different  in  both  cases/6*  If  the  material  can 
be  described  as  a  homogeneously  broadened  two  level  system  then  the  variation  of  the 
absorption  coefficient  with  intensity  is  given  by  : 


a((o)  = 


«n(a>) 

u-L 


(3.1) 


where  a0(a)),  /  and  Isat  are  the  unsaturated  absorption  coefficient,  the  input  intensity 
and  the  saturation  intensity  respectively. 


For  extremely  inhomogenously  broadened  systems  the  expression  for  the 
absorption  coefficient  is  modified  to  read  : 


a(o)) 


(3.2) 
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Ill  physical  lams,  ilia  lnhomogcnotiN  pain  01  absorption  salinates  iiiau1  slowly 
.it  large  intensities  because  the  signal  can  draw  energy  I'm  in  an  increasingly  w  uler  range 
of  packets  as  the  signal  intensity  increases. 


liquations  .VI  and  3.2  hold,  provided  that  die  contribution  from  the  real  part  of 
the  susceptibility  can  he  ignored,  i.e.,  there  should  he  no  aperture  present  in  front  of  ih  • 
signal  photodiode. 

The  experiments,  with  no  aperture,  were  performed  on  10  different  samples  and 
the  resultant  data  was  fitted  by  a  nonlinear  least  squares  lit  using  both  equations  3. 1  and 
3.2.  with  lM  being  the  only  fit  parameter.  In  all  the  eases,  equation  3.2  (i.e.  the 

equation  for  inhomogenously  broadened  systems  )  gave  a  far  superior  fit.  A  typical 
result  is  shown  in  figure  3.2.  The  sample  in  the  ligure  is  an  I8nm  colloid  dispersed  in 
deionised  water. 


Intensity  (\V/cm  ~) 

l’iguic  3.2  :  Absorption  Coefficient  Vs  Iaiser  intensity.  A  -  522/imi 

From  the  least  squares  fit  it  is  possible  to  determine  /J<J(  which  is  then  used  in  the 
equation 

Im  y(5>  = 

*  ~>nl 

*«  (3.3) 

to  calculate  the  imaginary  part  of  the  third  order  susceptibility. 

/i0  and  are  the  permittivity  of  free  space,  the  linear  index  of  refraction  and  the 
low  intensity  absorption  coefficient  of  the  sample  respectively. 
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liquation  3.3  can  ho  derived  by  Marline  lion)  Maxwell's  equations  as  lollou.s  : 


l;or  nonmagnetic,  centrosym metric  materials,  the  polarisation  induced  by  the  incidence 
of  an  intense  liehl  source  is  given  bv 


p  =  ^{x"^/:me.e)e 


l  Ising  this  expression  in  the  equation  for  the  displacement  vector 


D  =  i'E  +  P 


we  get 


D  =  £„(l  +  7"'  +  Xa'E.E)E  =  e0(e,+e„)E  (3  6) 

where  £,  =  1  +  %n)  and  £nl  =  x(})E.E,  are  the  linear  and  nonlinear  parts  of  the 
permittivity  of  the  medium. 

The  refractive  index  of  the  material,  n,  will  also  undergo  a  change  due  to  the  incident 
light.  This  is  expressed  by 


n-nQ  +  n2I 


where  n0  and  are  the  linear  and  nonlinear  indices  of  refraction. 


Refractive  index  is  related  to  the  permittivity  of  the  medium,  £,  through 


e  =  n2  =  +  2n</t2/  +higher  order  terms 

Here  n]  =  e,  and  In^l  =  £nl. 


Thus, 


2  nQn2l  =  X0)E.E 


which  implies  that. 


2  EqCIIq 


(3.10) 
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For  an  absorbing  maicrial  wo  have  the  equations 

/  =  /„  e.\p( -(c)  and  >'  ~  K, oxp (ink:.)  ,  V  ,  2 ) 

where  /„  and  /:„  are  the  incident  light  intensity  and  electric  field  respectively,  a  is  the 
linear  absorption  coefficient  and  k  is  the  wave  vector. 

From  equation  3.  1 2  we  get 


a  =  -lm(n)A:  = 


I  m  x  (l) 

2colh0  ( 


(3.13) 


The  term  a  is  the  slope  of  the  linear  pari  of  the  intensity  dependent  absorption  line.  It 
can  be  approximated  to 


(3.14) 

2/* 

Combining  equations  3.12  and  3.13  results  in 
__  *'c)cr,n  M 

2k/’«  (3.15) 

which  is  the  same  as  equation  3.3. 

Once  /,„  has  been  determined  from  the  nonlinear  least  squares  fit  and  an  from 
a  low  intensity,  linear  absorption  spectrum,  it  is  a  simple  matter  to  compute  the 
imaginary  pan  of  the  third  order  susceptibility.  This  was  done  for  all  of  the  samples  and 
the  results  are  tabulated  in  table  3. 1 . 

The  real  pan  of  the  susceptibility  can  be  measured  by  placing  an  aperture  in 
front  of  the  signal  detector,  as  described  earlier  on  in  this  chapter.  It  was  found  that  the 
trace  observed  in  this  case  did  not  exhibit  increased  transmission  peaks  and  decreased 
transmission  valleys,  indicative  of  a  refractive  index  change  but  was,  in  fact,  identical 
to  the  trace  obtained  with  no  aperture.  This  implies  that  at  these  wavelengths,  that  are 
very  close  to  the  surface  plasmon  resonance  frequency,  it  is  the  imaginary  part  of  the 
susceptibility  that  dominates  over  the  real  part.  This,  of  course,  comes  as  no  surprise 
because  ihc  real  part  of  (he  susceptibility  is  expected  to  be  low  in  the  vicinity  of  a 
resonance  decreasing  to  zero  at  its  peak.  The  imaginary  part  of  the  susceptibility. 
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should,  on  the  other  hand  increase  to  its  maximum  value  at  the  peak  ol  the  resonance 
frequency. 


RESULTS  OBTAINED 


Sample 

Size 

A 

a0(an~') 

516nm 

«o  («»"') 

522nm 

‘s, 

W/ am 

516nm 

— 

Is, 

Wirin' 

522nm 

nrV': 

516nm 

m'V" 

522nm 

N 

395 

5.50 

5.74 

3.45x  10s 

3.45x10s 

1 . 1 2x  Hf19 

6.49x1  (f20 

M 

392 

8.81 

9.32 

4.70x10s 

1.89x10s 

7.23x10'’° 

1.921 0'19 

() 

365 

6.10 

6.29 

1.05x10s 

1.10x10s 

2.24xl0'w 

2.23.x  10'19 

R 

326 

25.79 

26.82 

3.18x10s 

2.36x10s 

3.13xl0'19 

4.43x  10"19 

S 

184 

4.67 

4.78 

8.01xl07 

3.75x10s 

2.25x  10‘2° 

4.97xl0'20 

Y 

179 

9.84 

9.91 

1.93x10s 

1.52x10s 

1.97xl0'19 

2.54xl0'19 

mm 

172 

15.94 

15.98 

2.72x10s 

1.86x10s 

2.26x  10~19 

3.35xl0'19 

z 

164 

4.69 

4.75 

1.63x10s 

1.81x10s 

1.1  lx  10-19 

1.02x  10"19 

p 

52 

8.69 

8.82 

5.32x10s 

7.13x10s 

6.3xlO~20 

4.82x  10-2° 

Q 

47 

8.81 

9.09 

2.78x10s 

1.40x10s 

1.22xl0'19 

2.53x  10"19 

Table  3.1  :  Experimental  Results 


By  observation  of  the  table  it  is  obvious  that  there  is  no  relationship  between  the 
size  of  the  inclusions  and  the  nonlinearity  of  the  solution.  What  is  evident  is  that 
Im£t3)  is  proportional  to  the  square  of  the  absorption  coefficient.  This  result,  at  a 
wavelength  of  522nm,  is  shown  in  figure  3.3. 
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0  5  10  15  20  25  30 

Absorption  Coefficient  a  tan  1 ) 

Figure  3.3  :  lin*(3)  Vs  a„.  X  =522 nm. 

At  higher  concentration,  the  value  for  lm*i3)  seems  to  level  off.  Past  a  certain 
volume  fraction  of  particles,  the  Maxwell  Garnett  approximation  for  calculating  the 
local  field  in  the  inclusions  is  invalidated.  If  the  particles  are  sufficiently  close  to  interact 
significantly  with  each  other,  the  field  in  each  individual  one  is  reduced,  leading  to  a 
correspondingly  lower  value  for  the  macroscopic  susceptibility. 

The  quantity,  Im^<3),  in  the  above  table  and  graph,  refers  to  the  nonlinearity  of 
the  solution,  not  that  of  the  actual  gold  panicles.  In  order  to  extrapolate,  from  this,  the 
correct  value  lor  the  inclusions,  the  concentration  of  the  solutions  and  the  local  field 
(actor  discussed  in  chapter  2  have  to  be  taken  into  account. 

l:rom  equation  2.42  it  can  be  seen  that  the  susceptibility  for  die  inclusions,  *,l3\ 
is  related  to  the  macroscopic  susceptibility  of  the  solution  through  : 

yO) 

vO>  _  X 

"H/,17?  '  6) 

where  p  is  the  volume  fraction  of  the  inclusions  and  /,  is  the  local  field  factor.  Using 
this  equation  it  is  possible  to  calculate  X?]  and  results  are  tabulated  in  table  3.2. 


4X 


Sample 

(A) 

Im*1" 

5 1 6nm 

,nV~2 

*in  X 
522nm 

nrV'2 

P 

xl()  - 

/. 

516nm 

/, 

522nm 

516nm 

m'V~2 

522nm 

nrV 2 

N(395) 

-1.12x10““* 

-6.49x10"" 

1.90 

0.992 

1 .090 

6.06x  ur15 

2.4lxl(rl:’ 

M(392) 

-7.23x10  u 

- 1 ,92x  10~‘9 

3.16 

2.65x  10~15 

4.81x  10-li 

0(365) 

-2.24xl0’19 

-2.23x  1()'19 

BW 

0.937 

1.031 

1.31x1 0-14 

8.88x10’" 

R(326) 

-3.13xlO"19 

-4.43x  10'19 

9.98 

0.886 

0.975 

5.09x  1()-15 

4.91x  1()’15 

S(1 84) 

-2.25x  1()’20 

-4.97x1  O'20 

1.14 

1.114 

1 .222 

1 .28x  l()“,s 

1.95x10’" 

Y(  179) 

-1.97x  10’19 

-2.54x  10’19 

2.65 

1.257 

1.378 

2.97x  1()“15 

2.66x10’" 

K(172) 

-2.26x  10’19 

-3.35x  10'19 

4.37 

1.230 

1.350 

2.26x  1()"15 

2.11x10’" 

Z(164) 

-1.1  lxlO’19 

-1.02x  10'19 

1.28 

1.247 

1.367 

3.59x  10~15 

2.29x10’" 

P(52) 

-6.30x  10’2° 

-4.82x  10’2° 

3.13 

0.942 

1.036 

2.56x  10"15 

1.34x10’" 

Q(  47) 

-1.22xl0-19 

-2.53xl0’19 

2.77 

1.085 

1.191 

3. 17x  10”15 

4.53x10’" 

Table  3.2  :  Values  of  x?]  calculated  from  the  macroscopic  susceptibility  of  the 
solution. 


CONCLUSIONS 

The  most  important  conclusion  to  be  drawn  from  these  results  is  that  xl3)  is 

also  independent  of  the  size  of  the  particle,  at  least  within  the  size  range  of  this  study. 
The  values  obtained  for  are  almost  an  order  of  magnitude  higher  than  those 

obtained  by  Hachc  et  alP-  8-  This  is  not  surprising  since  they  calculated  their  local 
field  factor  by  comparing  the  absorption  spectrum  of  their  largest  particle,  which  they 
approximated  to  have  bulk  values,  to  the  absorption  spectrum  of  all  the  other  particles, 
to  obtain  e, ,  while  in  this  study,  die  absorption  spectrum  of  the  particles  was  compared 
to  an  actual  bulk  gold  spectrum,  to  calculate  t',  The  latter  method  yields  a  lower 

value  for  the  local  field  factor  and  hence  a  higher  valve  for  the  intrinsic  susceptibility  of 
the  inclusions.  By  way  of  comparison,  it  is  interesting  to  note  that  Bil,  microclustcr 
colloidal  solutions  are  reported  to  have  £,l3)  values  almost  two  orders  of  magnitude 
greater  than  these  ones  for  gold  clusters/6) 

If,  within  experimental  error  £,0)  can  be  taken  to  be  constant,  it  turns  out  that 
the  susceptibility  of  the  solution  is  proportional  to  the  fourth  power  of  the  local  field 
factor  which  is  in  turn  proportional  to  the  square  of  the  absorption  coefficient.  This  is 
exactly  the  same  as  the  result  expressed  in  figure  3.3. 

In  conclusion,  using  a  simple  but  effective  single  beam  method,  it  is  possible  to 
evaluate  the  sign  and  magnitude  of  both  the  real  and  imaginary  pans  of  the  nonlinear 
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nomogcncousl  v  01  inhomogenousiy  broadened.  Ilk'  mhomogcnousl  v  broadened 
colloidal  gold  solutions  used  in  litis  study  showed  a  nonlinear  response  dominated  hv 
the  imaginary  part  of  the  susceptibility,  which  is  to  expected  since  the  experiments  were 
jviiormed  in  the  vicinitv  ol  the  surface  plasmon  resonance  frequence  Die  ellective  or 
macroscopic  nonlinearity  ol  the  solution  is  negative.  This  is  implied  by  the  fact  that  the 
absorption  decreases  with  increasing  intensity.  The  nonlinearity  ol  the  inclusions  is 
related  to  the  nonlinearity  of  the  solution  through  the  local  field  factor  l\ .  At  these 
wavelengths.  J  ~  is  nearly  real  and  negative  which  leads  to  the  conclusion  that  jr;'1  is 
nearly  imaginary  and  positive.  Thus,  while  the  absorption  ol  the  colloidal  solution 
decreases  with  intensity,  the  absorption  of  the  inclusions  actually  increases!  This  is  a 
consequence  of  the  local  field  correction. 

Despite  the  fact  that  saturated  absorption  spectroscopy  is  a  very  powerful  and 
useful  technique,  it  provides  no  information  on  the  dynamics  of  the  interaction  between 
the  nonlinear  material  and  the  electromagnetic  light  field.  To  do  this  at  least  two  or  more 
beams  are  required.  The  following  three  chapters  describe  different  experimental 
techniques  in  which  this  is  done. 
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CHAPTER  4 


DEGENERATE  FOUR  WAVE  MIXfNfi  USING  A  COHERENT 

LIGHT  SOURCE 


Four  wave  mixing  is  one  of  the  most  interesting  and  thoroughly  investigated 
nonlinear  optical  effects/ It  is  a  process  that  is  allowed  in  all  media  and  has  therefore 
found  numerous  important  applications.  Because  of  the  inherent  resonant  enhancement 
feature,  it  has  attracted  a  great  deal  of  attention  as  a  modem  spectroscopic  technique. 
With  appropriate  arrangements  four  wave  mixing  can  be  used  to  study  transitions 
between  excited  states,  to  measure  longitudinal  and  transverse  relaxation  time  etc. 
Thus,  experiments  on  transient  four  wave  mixing  can  be  used  to  study  not  only  the 
particular  time-ordered  four  wave  mixing  process  itself,  but  also  to  obtain  the  various 
relaxation  rates  in  the  medium. 

In  general,  transient  four  wave  mixing  deals  with  the  situation  where  three  input 
pulsed  fields,  of  either  the  same  or  different  frequencies,  interact  in  a  medium  in  a  given 
time  order.  The  radiative  output,  which  is  the  fourth  wave,  shows  a  time  variation 
depending  on  the  time  sequence  and  separation  of  the  input  pulses.  If  the  four  waves 
are  all  at  the  same  frequency,  co,  the  interaction  is  termed  degenerate  four  wave 
mixing,  hereafter  called  DFWM.  In  general,  two  of  the  fields  arc  strong  and  are  called 
the  pump  fields  and  the  weaker  third  field  is  called  the  probe  field. 

The  second  nonlinear  polarisation  term  in  equation  1.1,  which  expresses  the 
relationship  between  polarisation  and  the  applied  field, 

P  =  £0(XiUE  +  Xa)E2  +  £(3,E3+ . \  (4l) 

i.e.,  P  =  £0Z(3) E.E.E  is  the  one  responsible  for  this  process.  (Recall  that  for 
centrosymmetric  or  non-isotropic  materials  xi2)  =  0). 

If  the  two  pump  beams  are  counter  propagating  and  the  probe  beam  enters  at  a 
small  angle  with  respect  to  them,  as  shown  in  figure  4.1,  the  signal  beam  propagates 
backwards  through  the  medium  and  remains  everywhere  a  phase  conjugate  replica  of 
the  probe  beam.  This  process  is  called  "Phase  Conjugation".  It  is  technologically  very 
important  because  it  is  possible  to  use  it  to  design  self-adaptive  optical  systems  that 
compensate  for  time  varying  phase  distortions  in,  for  example,  high  gain  laser 
oscillators  and  amplifiers,  optical  fibres,  etc/2) 
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I  *11111 [>  2 

•+ - 

Pmbe 

,> '  Signal 

figure  4. 1  :  Sdicmaiic  diagram  ol  a  phase  conjugation  experimental  .set-up 

The  first  observations  ol  phase  conjugation  bv  lour  wave  mixing  were  reported 
by  Bloom  and  Bjorkland*^  and  Jensen  and  llellwarth*4),  both  parties  using  CS?  as  the 
nonlinear  medium. 

THEORETICAL  FORMALISM 

The  fields  involved  in  the  interaction  can  be  described  as  plane  waves  using^); 


Bump 


Sample 


E,(r.t)  =  j{AM, )exp|i(*,r  -  n>/)|  +  c.c.}  (4.2) 

where  k{  is  the  complex  propagation  vector,  rj  is  the  distance  along  k,  and  c.c.  denotes 
the  complex  conjugate.  A-,  can  be  described  as 

kt  =  *,(*'  +  ia  /  2)  (4.3) 

where  kt  is  a  unit  vector  in  the  direction  of  propagation  of  a  is  the  intensity 
absorption  coefficient  in  die  medium  at  the  f  requency  r o  and  k:  = 


Thc  induced  nonlinear  polarisation  can  similarly  be  wntten  as  : 
PNl'(r,t)  =  ^{P(r)exp[-/<y/J  +  c.c.} 


(4.4) 


Where, 


p(r),  =  Tfc'oX{z^(-"-W,£U,-w):/\,(r|)  Aj(rJ  A;(r,)  cxpfii-.f*,  +  L,  - k'3 

4  iU  K  v 


(4.5) 


By  substituting  for  PSI  in  the  inhomogenous  wave  equation. 
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<)'!■:  ,a>s'  ,  t 

p-  <4<>) 

<?/•  (Jr 

and  applying  the  slowly  varying  envelope  approximation  (SVEA)  i.c., 

k2At  »  T-  (4.7) 

<7C  f  A: ' 

it  is  possible  to  solve  lor  the  field  amplitudes  Aa  and  A\  which  represent  the  new 
"signal"  beam  and  the  probe  beam  respectively.  By  standard  techniques,  the  following 
equations  are  obtained. 


=  iQAl  cxp(-ccz) 
dz 


=  iQAtcxp{az) 


The  coupling  constant  Q  is  given  by 


~Xl3)A,  (0)  Aj  ( £)expj^-  ^  aL 


(4.10) 


where  n  is  the  complex  index  of  refraction  and  L  =  tt— rv,  where  L  is  the  sample 

M 

thickness. 

Using  the  boundary  conditions  Aa(L)  =  0  (  no  phase  conjugate  wave  initially  )  and 
4(0)  =  A 3  ( the  initial  value  of  the  probe  Held  ),  the  solutions  of  equations  4.8  and 
4.9  are(5\ 


4(Z)  = 


2 iQA] sin[//(Z  -  L) /  2]exp(-oZ /  2) 
asin(HL  /  2)  +  Hcos(HL  /  2) 


(4.11) 


4(z)  = 


-4{asin[//(Z  -  L) /  2]  -  //cos[//(Z  -  L) /  2]}exp(aZ /  2) 


as\n(HL  1 2)  +  Hcos(HL  /  2) 


(4.12) 


//  =  [4|<2J2 


(4.13) 
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I  Ik'  phase  conjugate  ivllccuv  1 1 s  is  delmed  as  the  ratio  ot  t Ik*  micnsiu  ol  the  generated 
signal  Ivam  lo  that  ol  the  prolre  Ircam.  It  i.s  given  hy. 


.4,(0) 

2(Tsin(  ///.  /  2) 

-<(<>) 

irsinl ///.  /  2)  +  //cos(  HU  2) 

In  the  limit  ol  low  reflectivity. 


R  =  ^i-[l  -exp( -<//.)]' 


(4.15) 


The  third  order  .susceptibility,  x  ls  related  to  the  relleetivity,  R  through'-*'- ()-  7) ; 


ot  _  4cVc0aV7e 

X  3fl)/Vf(l-T) 


(4.16) 


where  i.s  a  the  linear  absorption  of  the  sample,  and  T  i.s  the  transmission  through  it. 


EXPERIMENT 

Phase  conjugation  is  only  one  of  the  many  different  geometrical  configurations 
for  DFWM.  Another  arrangement  that  is  commonly  used  is  called  the  "folded  boxcar" 
configuration. 9)  It  was  the  one  used  in  this  study  and  is  shown  in  figure  4.2.  It 
consists  of  two  pump  beams  that  overlap  and  interfere  in  the  sample  to  produce  an 
intensity  modulation  that  is  subsequently  detected  by  a  probe  beam.  This  intensity 
modulation  can  affect  the  sample  ;n  many  different  ways  :  it  can  result  in  a  refractive 
index  modulation,  a  population  lodulation.  absorptive  index  modulation  etc.  The 
modulation  can  be  visualised  as  a  diffraction  grating  within  the  sample,  i.e..  the  pump 
beams  "write”  a  grating.  "Reading"  this  grating  can  give  information  on  different 
physical  properties  of  the  medium.  This  i.s  done  by  the  weak  probe  beam.  When  the 
probe  beam  is  incident  on  the  grating  it  diffracts  in  a  new  direction  giving  rise  to  the 
signal  beam.  The  signal  beam  i.s  detected  by  a  calibrated  photodiode/lock-in-amplificr 
(LIA)  combination.  The  L1A  is  triggered  by  a  mechanical  chopper  in  one  of  the  incident 
beam  paths. 
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Figure  4.2  :  DFWM  using  the  "folded  boxcar"  configuration.  (C.D.  :  Cavity 
Dumper) 


It  is  now  usual  to  talk  about  the  diffraction  efficiency  which  is  defined  in  the  same 
manner  as  the  phase  conjugate  reflectivity. 

The  relationship  between  x°]  and  the  diffraction  efficiency  rj  is  given  by 


m  =  8  c2n2e0ajn 
X  3ffl/Vf(l  -  T) 


(4.17) 


Thus,  in  order  to  determine  the  third  order  susceptibility  of  a  particular  sample 
one  need  only  measure  the  ratio  of  the  energy  of  the  diffracted  beam  to  that  of  the  probe 
beam.  It  is  possible  to  separate  out  the  different  tensorial  components  of  x0)  by 
varying  the  polarisations  for  the  pump  and  probe  beams.  By  introducing  a  time  delay 
between  the  probe  beam  and  either  of  the  two  pump  beams  it  is  possible  to  determine 
the  relaxation  time  associated  with  the  grating  and  hence  the  decay  time  of  the  physical 
property  that  was  modulated  by  the  interference  pattern  of  the  pump  beams. 

The  laser  system  used  in  this  study  was  a  cw  mode  locked  Ndu  doped  YAG 
laser  operating  at  a  repetition  rate  of  100MHz.  The  output  of  the  laser  was  frequency 
doubled  and  then  used  to  synchronously  pump  a  dye  laser.  The  dye  used  was 
Rhodaminc  6G  which  provided  a  tunability  from  560  -  61()nm.  The  output  from  the 
dvc  laser  was  fed  into  a  cavity  dumper  in  order  to  reduce  the  repetition  rate.  It  was 
possible  to  vary  the  repetition  rate  from  100MHz  down  to  about  40KH/..  The  width  of 
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1 1  ic  d\o  laser  pulse  was  mcusuied  In  autocorrelation  ieehnu|ues  iisin::  a  second 
harmonic  generating  crystal  and  found  to  (>e  ~5ps.  A  cross  correlation  of  the  two  pump 
beams  was  taken  at  the  sample  to  ensure  that  they  arrived  at  that  point  simultaneously. 
The  delay  of  the  probe  beam  was  then  scanned  to  find  which  setting  produced  the 
maximum  signal. 

liven  at  the  lowest  repetition  rate  of  40  KH/.  the  only  sample  that  it  was  possible 
to  do  measurements  on  was  a  5()A  gold  colloid  dissolved  in  acetone. 

The  sample  was  examined  at  a  range  of  different  repetition  rates  at  a  wavelength 
of  58()nm.  The  concentration  was  the  same  in  all  cases.  The  results  are  plotted  in  figure 
4.3.  The  background  at  r  <  0  is  just  scatter  from  die  chopped  beam  that  is  collected  bv 
the  signal  detector.  This  is  verified  by  blocking  the  other  two  incident  beams.  From  the 
figure  it  can  be  seen  that  as  the  repetition  rate  is  increased,  the  signal  decays  with  a 
continuously  slower  relaxation  time.  This  is  indicative  of  a  thermal  build-up  in  the 
solvent  due  to  absorption  of  energy  by  the  solute.  For  the  water  soluble  colloids,  the 
thermal  build  is  more  drastic  making  it  impossible  to  perform  any  meaningful  time 
resolved  experiments,  even  at  the  lowest  repetition  rate  of  the  laser. 

When  all  the  three  input  beams  ate  lime  coincident  at  the  sample,  the  magnitude 
of  the  DFWM  signal  is  related  to  |^<3)|  .  It  can  be  determined  by  measuring  the 

diffraction  efficiency  of  the  sample,  relative  to  that  of  a  well  characterised  solvent.  This 
is  done  to  allow  for  variations  in  the  beam  quality.  The  solvent  in  this  case  was  CS, 
which  has  a  x ^  value  of  b.8x  l()'|,esu.<l0)  The  tensorial  component  refers  to 
the  case  when  all  three  beams  are  linearly  polarised  in  the  same  direction,  which 
corresponds  to  the  experimental  situation. 

When  the  measurements  are  done  relative  to  a  reference  sample,  x"'  *or  losl 
material  is  obtained  using  the  formula. 


vO)  _  v(3) 

a  -  Xcs,  ■ 


/ 


A-., 


'os, 


acxp(cui  /  2) 
1  -cxp(-cui) 


(4.18) 


where,  /  and  are  the  diffraction  efficiencies  and  n  and  nCSj  are  the  linear  refractive 
indices  for  the  sample  and  for  the  reference  material,  CS2 ,  respectively,  and  a  is  the 
linear  absorption  coefficient  of  the  sample. 
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Delay  Tijne  (ps) 


RESULTS  OBTAINED 

It  is  possible  to  examine  the  relationship  between  the  absorption  coefficient  and 
susceptibility  by  performing  the  experiment  at  a  range  ol  different  wavelengths.  The 
wavelength  was  varied  from  562-606nm.  and  the  results  obtained  are  shown  m  fable 
4.1. 


Wnvi  lt'ti^ili  (nm) 

A  bsorption  Coefficii  nr 
(cm  ') 

C'K'i'') 

562 

50.19 

1.62x10''" 

560 

46.04 

1. 29x10'“' 

575 

42.72 

1.00x1  O'19 

581 

39.28 

8.67x  10'2° 

587 

36.48 

6.95xl0'20 

594 

32.95 

4.58x  10"’° 

600 

30.46 

3.53xlO~20 

606 

28.21 

3.24xl0'20 

Tabic  4.1:  Results  obtained  for  wavelength  dependent  measurements 


The  results  are  shown  graphically  in  figure  4.4. 


2.5  3  3.5  4  4.5  5  5.5 

Absorption  Coefficient  (mm1) 
f  igure  4.4  :  x'"  Vs  Absorption  Coefficient,  a 


As  expected,  x  ’  is  again  proportional  to  the  square  of  the  absorption  coefficient.  The 
tailing  off  of  x  w'ih  increasing  absorption  coefficient,  as  observed  in  chapter  3.  is 
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not  scon  here  as  the  concentration  is  well  below  the  threshold  value  lor  the  Maxwell 
Garnett  approximation,  lor  the  local  field  in  the  inclusions,  to  be  applicable. 
Furthermore,  the  absorption  is  changed,  not  by  varying  the  concentration,  but  by 
varying  the  wavelength. 

The  value  for  the  susceptibility  of  the  inclusions  is  obtained  from  the 
macroscopic  susceptibility  using  the  formula  2.42, 


3)  _  _ X 


(3) 


xl3)  = 


M  ft 


(4.19) 


The  results  are  tabulated  in  table  4.2. 


/. 

1.62xl0"19 

1.183 

4.02xl()"16 

1.29xl0~19 

1.138 

3.74x  10~16 

l.OOxKT19 

1.092 

3.42x1 0"16 

8.67x  10-2° 

1.040 

3.60x  10-16 

6.95xlO-20 

0.990 

3.52xl0“16 

4.58x  10-20 

0.928 

3.00xl0~16 

3.53X10-20 

0.877 

2.90x  ur16 

3.24x  10'2° 

0.828 

3.35x  1()'10 

Table  4.2  :  values  for  the  susceptibility  of  the  inclusions,  taking  into  account  the 
local  field  factor,  /, .  The  value  for  p  is  2. lx  10"3 . 


CONCLUSIONS 


These  results  imply  that  the  susceptibility  of  the  inclusions  is  independent  of  the 
wavelength  or  the  absorption  coefficient  of  the  colloidal  solution.  The  only  quantity  that 
varies  with  wavelength  is  the  local  field  factor  /,.  From  equation  2,y  it  can  be  seen  that 
the  absorption  coefficient,  a,  is  directly  proportional  to  the  square  of  the  local  field 
factor.  Since  xa)  is  proportional  to  the  local  field  factor  to  the  fourth  power,  ( taking  p 
and  x{3)  to  be  constant,  within  experimental  error  ),  it  should  also  be  proportional  to 

the  square  of  the  absorption  coefficient.  This  is  the  result  shown  in  figure  4.4.  The 
values  obtained  for  £,U)  in  the  wavelength  range  562-606nm  are  of  the  same  order  of 

magnitude  as  those  obtained  by  Hache  et  al.  at  resonance/11-  12-  13>  The  reason  for 
this,  as  explained  in  chapter  3,  lies  in  the  value  of  the  local  field  factor  and  the  method 
that  is  used  to  determine  it. 
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I'Ik*  \  .linos  in  lik"  above  table  .no  absolute  ■cnee.  m  ihi.s  expei imemai 
arrangement.  there  is  no  wav  of  determining  the  sign  »’i  /  .  Ii  is  also  not  possible  lo 

say  what  contribution  ol  the  diffracted  signal  comes  from  the  real  part  of  the 
suseeptihilitv  and  what  comrihulion  comes  from  the  imaginary  part.  The  local  field 
factor  I,  at  these  wavelengths  has  real  and  imaginary  parts  o!  equal  magnitude.  both 
being  positive.  This  information  is,  however,  redundant  since  the  magnitudes  of  the 
real  and  imaginary  parts  of  the  macroscopic  susceptibility  are  unknown.  This  means 
that  the  relative  magnitudes  of  the  real  and  imaginary  parts  olThe  susceptibility  of  the 
inclusions  also  remains  unknown.  The  only  thing  that  c<;/i  Iv  said  is  that  the  real  part  of 
the  susceptibility  is  probably  not  negligible.  1'his  was  the  case  in  Chapter  3  (  i.o..  the 
real  part  of  the  susceptibility  was  negligible  iwhere  the  experiments  were  being 
performed  very  close  to  the  surface  plasmon  resonance  frequency.  As  explained  in 
Chapter  3,  the  further  away  from  resonance  that  the  experiments  are  performed,  the 
greater  should  be  the  relative  contribution  of  the  real  pan  of  the  susceptibility. 

I'he  change  in  the  DFWM  signal  as  a  function  ol  time  delay  of  the  probe  beam 
gives  the  decay  of  the  grating  and  hence  the  time  response  of  the  optical  nonlinearity. 
The  experiments  were  all  performed  at  a  repetition  rate  of  286  H/..  The  response  time  in 
cadi  case,  which  is  a  higher  order  autocorrelation  of  the  laser  pulses  convoluted  with 
the  material  response,  was  found  to  be  shorter  than  the  intensity  autocorrelation  width, 
indicating  that  the  experimental  resolution  is  limited  by  the  laser  pulse  width.  The  only 
information  to  be  gleaned  from  this  experiment,  therefore,  is  that  the  lifetime  of  the 
excited  state  in  the  resonant  two  level  system  is  ai  most  5ps.  This  result  is  consistent 
with  measurements  performed  by  Hache  <•;  a/,*11-  Bloomer  <7  a/*14)  and  Heilweil 
and  Hochstrasser.C-*'!  A  similar  result  was  found  by  yang  ei  ol  for  Copper  nanoclustcr 
colloidal  solutions.* 

The  relaxation  processes  in  a  resonant  two  level  system  arc  usually  described  by 
a  second  parameter,  7\  along  with  the  lifetime  of  the  excited  state.  Tr  also  known  as 
die  longitudinal  relaxation  time.  The  parameter  T:  is  known  as  the  transverse  relaxation 
time  or  the  phase  relaxation  time  .  It  represents  the  decay  of  the  coherence  between  the 
two  levels.  By  an  appropriate  arrangement  it  is  possible  to  determine  T:  from  the  above 
DFWM  experiment  as  well.  However,  T2  generally  lies  in  the  fscc  time  regime  and  so, 
using  coherent  5  psec  pulses  would  give  very  little  information  about  it.  The  following 
chapter  describes  a  novel  technique  for  determining  T:  using  incoherent  light. 
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CHAPTER  5 


DEGENERATE  FOUR  WAVE  MIXING  USING  INCOHERENT 

LIGHT 


In  order  to  understand  the  dynamical  behaviour  of  light-matter  interaction,  it  is 
very  important  to  study  the  relaxation  processes  associated  with  the  excited  stales  of 
materials.  In  many  cases  of  interest,  especially  in  condensed  matter,  these  relaxation 
processes  take  place  on  a  subpicosecond  or  femtosecond  timescale.  In  transient 
spectroscopy  using  coherent  light,  the  time  resolution  of  the  experiment  is  limited  by 
the  pulse  width  of  the  exciting  radiation.*1  • 21  Thus  in  order  to  probe  the  dynamical 
behaviour  of  such  materials  effectively  it  is  essential  to  have  access  to  ultrafast  lasers. 
To  this  end,  a  lot  of  effort  has  been  devoted  to  obtaining  short  pulse  lasers  and  with 
CPM  dye  lasers,  pulses  with  widths  of  less  than  100  fs  have  been  obtained.  Even 
shorter  pulses,  with  pulse  widths  of  6  fs  have  been  generated  by  using  a  pulse 
compression  technique.*3’  41  There  are,  however,  several  problems  associated  with 
using  such  short  pulses.  Firstly,  the  optics  required  to  generate  or  compress  ultra  short 
pulses  are  complicated  and  prohibitively  expensive.  Also,  as  there  are  not  very  many 
suitable  combinations  of  laser  dye/saturable  absorber,  ultra  short  pulses  can  only  be 
generated  in  a  very  limited  wavelength  regime.  Finally,  even  if  one  were  able  to 
succeed  in  getting  a  short  enough  pulse  at  a  suitable  wavelength,  there  arc  further 
problems  associated  with  the  dispersion  of  the  pulse  as  it  traverses  through  the  optics  of 
the  experiments. 

One  solution  to  this  problem  is  to  perform  the  experiments  in  the  frequency 
domain.  In  1978  Yajima  et  al  *5-  proposed  resonant  Rayleigh-type  mixing 
spectroscopy  for  the  measurement  of  the  longitudinal  relaxation  time,  7\,  and  the 
homogenous  transverse  relaxation  time,  Tv  and  the  relaxation  times  in  some  materials 
have  actually  been  measured.*7)  However  in  this  case  the  analysis  of  the  results  is 
extremely  complex  and  can  lead  to  ambiguities  in  the  determination  of  these  limes. 

Prompted  by  these  difficulties,  in  1984,  Morita  et  a/*8)  and  Beach  and 
Hartmann  independently  came  up  with  a  new  spectroscopic  technique  that  utilises 
temporally  incoherent  light  to  determine  ultrafast  relaxation  times.  Traditionally  optica] 
coherent  experiments  have  always  been  performed  with  coherent  laser  light.  While  this 
may  he  essential  for  experiments  such  as  self  induced  transparency,  it  is  not  a  universal 
requirement.  Indeed,  not  only  can  optical  coherent  transients  be  generated  with  an 
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mcohorcm  light  .source,  iliciv  arc  decided  advantages  in  dome  so.  Optical  transitions 
can  only  he  induced  within  the  bandwidth  ol  the  excitation  source.  1  or  a  laser 
oscillating  in.  at  most,  a  lew  weil  defined  modes,  this  means  that  only  a  small  traction 
ol  the  atoms  or  molecules  within  the  resonance  line  participate  in  an  experiment.  This 
problem  can  he  overcome  by  performing  the  experiment  using  a  broadband  incoherent 
I mht  source  which  also  provides  ullrahigh  time  resolution.  The  basis  ol  this  method  is 
that  the  signal  light  generated  or  modulated  by  nonlinear  optical  effects  consists  of  a 
correlation  ol  the  material  response  and  the  input  field  and.  as  such,  is  dependent  on  the 
correlation  time  of  the  pulse  instead  of  its  width.  Temporally  incoherent  light  has  a  wide 
spectral  <  idlli.  At),  and  a  correspondingly  short  correlation  time  determined  by  I /At), 
which  is  much  shorter  than  its  temporal  width.  This  kind  of  light  appears  as  a  single 
pulse  of  duration  tc  in  an  autocorrelation  experiment ( l("  and  is  therefore  expected  to 

play  the  same  role  as  a  short  pulse  in  experiments  utilising  the  correlation  technique. 
This  fact  has  been  verified  experimentally  in  several  situations,  including  Raman 
spectroscopy/ 11  •  for  measuring  vibrational  relaxation  times,  and  Degenerate 

lour  wave  mixing  (DFWM)  (1 3- Id)  photon  echoes  bf  20)  and  pump  and  probe 
spectroscopy  (21). for  measuring  phase  relaxation  times  and  excited  state  lifetimes  of 
nonlinear  materials.  This  chapter  focuses  on  the  measurement  of  the  dephasing  time  of 
colloidal  gold  solutions  by  DFWM  spectroscopy.  In  the  self  diffraction  configuration, 

—4  -4 

two  temporally  incoherent  light  beams  with  wave  vectors  ki  and  lo,  originating  from  a 
single  beam  at  frequency  cu  are  made  to  overlap  in  a  nonlinear  sample/2)  Due  to  the 
third  order  nonlinearity  of  the  material,  two  new  output  beams  arc  generated  at  the  same 
frequency  but  in  the  directions  k 3  =  2 k2  -kt  and  k.  =  2k}  -  k,.  Even  if  the  pulse  width 
of  the  light  is  much  longer  than  either  7,  or  7\,  the  correlation  trace,  i.e.,  the  diffracted 
light  intensity  as  a  function  of  delay  ume.  r.  between  the  two  beams,  decays 
exponentially  with  a  lime  dependent  on  /',.  provided  that  the  coherence  time  of  the 
pulse,  rr,  is  less  than  7\/^  This  fact  is  extremely  useful  because  it  is  far  easier  to 
prepare  an  incoherent  light  source  with  a  short  correlation  time  rf  than  to  produce  an 
ultrashort  pulse  with  the  same  duration,  especially  in  the  range  less  than  1  ps.  In  an 
extreme  case,  even  cw  light,  which  has  an  infinite  duration,  is  capable  of  providing  fs 
time  resolution  provided  it  has  the  adequate  spectral  width.  In  the  event  that  rr  is  not 

very  short,  the  resolution  of  the  experiment  is  limited  on  it. 


THEORETICAL  BACKGROUND 


The  electric  fields  of  the  incoherent  light  beams  involved  in  the  interaction  can 
he  described  as  ^ 

E(rj)  =  jexp(-/mr  +  ik . r |  +  c.c.  (5. 1 ) 


with. 


$(t)  =  e{t)R(t)  (5.2) 

where  k  is  the  propagation  vector,  n  is  a  unit  vector  in  the  direction  of  k  ,  c.c.  denotes 
the  complex  conjugate,  e(t)  is  a  normal  function  and  R(t)  is  a  complex  random 
function  representing  a  stochastic  stationary  Gaussian  process,  which  is  how  the 
incoherent  beam  is  represented. 


For  such  a  function. 

( R\t)R(t  +  z))  =  f(z ) 

(5.3) 

(R(t)R(t+z))  =  (Rm(t)R-(t+rj)  =  Q 

(5.4) 

(/?*(0)  =  (/?('))  =  o 

(5.5) 

where  /( r)  is  the  correlation  function. 

The  total  electric  field  incident  on  the  material  is  a  superposition  of  the  two 
— *  — »  — * 

fields  with  wave  vectors  ki  and  k?  where  the  field  with  wave  vector  k2  is  temporally 

delayed  by  z  relative  to  the  other.  The  total  Held  can  be  written  in  the  form 


E(r,t)  =  E]  +  £,  =  | <*j  t  +  r-— —  exp(-/a>(r  +  z)  +  ik].r j  + 


n.r 

v 


exp(-t 


l(Ot  + 


+  c.c. 


(5.6) 


The  relaxation  processes  in  a  resonant  two  level  system  arc  usually  described  by 
the  two  parameters,  the  longitudinal  relaxation  time,  7j,  representing  the  decay  of  the 
population  difference  and  the  transverse  relaxation  time.  /’,  ,  representing  the  decay  of 
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ilk*  coherence  between  ilu*  two  levels.  Il  ilieie  is  .1  distribution  o|  uansiiitiii  ncquciicie-, 
between  the  levels  ilk*  system  is  considered  to  bo  mhomogeiiouslv  broadened 


Ilk*  output  light  field  at  tin*  I  requeue  v.  a),  m  a  lour  wave  mixing  process  is 
proportional  to  die  dual  older  induced  polarisation.'"’  which  has  the  lorni. 

/” '  (/  ./)  -  /’  ’(/  ,/)e\p(  -not)  1  <  ..  .  (> 

with. 

')  =  tVj  dtojijjy^^rj.to^to,,)  l5.<S) 

0 

where  N  is  the  atomic  number  density,  the  subscripts  u  and  b  denote  lower  and  upper 
levels,  respectively,  is  the  electric  dipole  matrix  element  of  the  transition  between 
the  two  levels,  and  g(ai0)  is  the  distribution  function  of  the  transition  frequency.  o>n, 
characterising  the  inhomogenous  broadening,  p)':t u>„)  is  the  third  order,  off- 
diagonal  density  matrix  element.  It  contains  four  components  at  the  wave  vectors  A,. 
A:.  A,  =  2A\  -  A,, and  A4  =  2A,  -  A,.  Since  the  experiment  involves  detecting  and 
measuring  the  signal  in  the  direction  A,,  it  is  only  necessary  to  consider  the  A, 
component  of  p^J . 

Tliis  is  given  hv  die  equation,'*) 

Phj[k,)=  -Up'0'  y  exp(/A:,./:  +  i(Ot)  X 

Jjr,  )  «//,  }<//,{$(/,)-!;(';)*’('»  +  r)exp[-/(w.,  -«>)(/,  -/,  f  t,)\ 

+  r)^(r,)exp[-/(w0  -w)(/, 

.Yexp[-y,(/,  -/:)-y?(rr  +h  -f,)|  (5.9) 

where  y,  =  T“'.  y,  =  7""',  p[0)  is  the  thermal  equilibrium  value  of  the  population 

difference  of  the  two  levels,  /,=/-  — — -  is  the  reduced  time  assuming  that 

v 

",  =  =  "  and  pha  =  =  // . 

OK 


The  ouipui  light  intensity,  wiiidi  is  the  quantity  that  is  actually  experimentally 
measured  is  proportional  to  |/>l3>(7:,/)|  .  However,  because  of  the  random  nature  of 
F"(f,t ),  it  is  necessary  to  calculate  the  statistically  averaged  value,  •/(£,),  of  this 
quantity,  which  can  he  written  as,(X* 


The  function  l:{x)  determines  the  shape  of  the  correlation  prollle.  It  can  be 
calculated  when  a  reasonable  form  is  given  to  the  correlation  function  /( r).  f(t)  is 
usually  approximated  to  be  a  delta  function.  Though  not  absolutely  correct,  this  choice 
is  reasonable  since  the  rate  to  be  measured  is  much  greater  than  the  coherence  time  of 
the  pulse.  In  the  event  that  the  coherence  time  is  comparable  to  or  greater  titan  the 
quantity  being  measured,  the  analysis  is  unnecessary  since  the  only  information  that  can 
be  obtained  in  such  a  case  is,  a  value  for  the  upper  limit  of  that  quantity. 

The  explicit  form  of  F{r)  when  g(m0)  is  a  Gaussian  distribution  centred  at  mn 
with  a  width  of  8(0  is  shown  for  two  extreme  cases  in  what  follows 

(i)  Homogenous  broadening  case.  So)  =  0. 

(i.i)  r>() 


F(t)  =  c, 


u  + 


— - — ^^exp(-2x)  -4(1-  u)exp[-(2  +  n).r] 


+  — ^exp(-2ux)  (5.11) 

1  -  u 

(i.ii)  r  =  () 

F(0)  =  c^2  +  2n  +  -U2j  (5.12) 

(i.iii)  T<0 

F(r)  =  c,[n  +  2exp(2x)]  (5.13) 


(ii)  Hxtremcly  inhomogenous  broadening  case.  So)  — >  °° 


Il  is  instructive,  at  this  stage,  to  compare  the  signal  behaviour  tn  this  ease  to  that 
in  a  conventional,  short  pulse,  two  beam  self  diffraction  experiment,  hereafter  referred 
to  as  short  pulse  DFWM.  The  similarity  between  the  processes  lies  in  the  fact  that  both 
profiles  reflect  the  relaxation  processes  governed  dominantly  by  the  phase  relaxation 
lime  T:.  There  arc,  however,  three  main  differences. 

(i)  The  correlation  traces  in  this  experiment  do  not  necessarily  decay  as  single 
exponentials,  while  those  of  short  pulse  DFWM  experiments  always  decay  single 
exponentially  at  a  rate  proportional  to  — . 
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(ii)  llio  present  correlation  profiles  have  non  zero  values  at  r  =  ±^,  while  no 
background  exists  in  the  short  pulse  DFWM  case. 

<iii)  In  this  experiment  the  correlation  profiles  at  r<0  grow  up  with  increasing  r  in 
the  homogenous  broadening  case,  while  those  in  the  short  pulse  DFWM  case  are 
always  zero  at  r<0. 

From  (i)  it  can  be  seen  that  the  profile  from  such  tin  experiment,  thus,  does  not 
necessarily  give  us  one  of  the  relaxation  limes  Tx  and  7',  in  as  unambiguous  a  manner 
as  die  usual  self  diffraction  experiments.  If.  however,  /',  »  7\ ,  the  profiles  represent  a 
single  exponential  decay  with  the  rates  277*  and  47V1  lor  the  homogenous  and 
extremely  inhomogenous  cases,  respectively.  These  rates  are  consistent  with  those  in 
the  short  pulse  case  and  T2  can  be  uniquely  determined.  In  condensed  matter,  such  as 
singlet-singlet  transitions  of  dye  molecules,  interband  transitions  of  semiconductors 
etc.,  T2  is  generally  very  short  and  often  falls  far  below  lps.  In  these  transitions,  the 
population  relaxation  time  Tx  is  often  much  longer  than  T2.  Therefore,  in  these  cases 
the  present  method  is  a  powerful  tool  for  determining  T2. 

In  the  event  that  7',  is  not  much  longer  than  T2,  the  present  method  cannot 
determine  the  relaxation  time  T2  definitely.  However,  even  in  the  worst  case,  it  is 
possible  to  determine  Tt  within  the  error  factors  of  1.5  and  2.7  in  the  homogenous  and 
extremely  inhomogenous  broadening  cases,  respectively.*8* 

EXPERIMENTAL  SET-UP 


The  experimental  set-up,  for  measuring  T;.  is  shown  in  figure  5.1.  The  PR  A 
dye  laser  was  pumped  by  a  PRA  Nitrogen  laser  and  emitted  5(X)ps  pulses  with  energies 
of  30  -  40  pJ  at  a  repetition  rate  of  1  -  5  Hz.  The  dye  used  was  Coumarin  485. 
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[  Ik*  output  beam  u.is  divided  into  two  pails  umii::  a  50:50  Ivam  spinier.  A I  lor 
puss  ini:  the  two  Ivams  through  variable  and  fixed  delay  lines,  they  were  reeombmed  in 
the  sample  (thickness  0. 1  mm)  using  a  20 cm  lens.  The  delay  time  t  between  the  beams 
was  changed  by  a  stepping  motor. 

The  spectral  width  of  the  emitted  radiation  was  increased  bv  replacing  the 
diffraction  grating  in  the  dye  laser  oscillator  cavity  with  a  mirror.'17-  11,1  I'nder  these 
conditions  the  bandwidth  was  found  to  be  10.5  nm  with  a  centre  wavelength  of  522nm. 
This  corresponds  to  a  coherence  time  of  Xb.5fs. 

Figure  5.2  shows  the  results  of  the  bandwidth  measurement.  The  experiment 
was  performed  by  passing  the  beam  through  a  monochromator  and  monitoring  us 
energy  at  the  exit  slit  as  a  function  of  die  wavelength  selected. 
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Figure  5.2  :  ( )utput  intensity  V\  A  .lor  measuring  Hie  bandwidth  ol  die  be;un 
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Hie  samples  studied  were  : 


RESULTS  OBTAINED  AND  CONCLUSIONS  ON  THE 
MEASUREMENT  OF  THE  DEPHASING  TIME 

The  correlation  profiles  lor  three  of  the  samples,  C,  S  and  M  are  shown  in 
figures  5.3,  5.4  and  5.5,  respectively.  These  three  samples  reflect  the  range  of  sizes 
used  in  this  study. 
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l;ipurc  5.3  :  Sample  Diffracted  intensity  Vs  delay  lime. 


73 
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figure  5.5  :  Sample  M,  Diffracted  intensity  Vs  delay  time. 

The  experiment  was  also  performed  on  a  dye  solution  of  Rhodamine  13,  as  a 
reference  material,  known  to  have  a  suh-ps  phase  relaxation  time/231  The  result  is 
shown  in  figure  5.6. 
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l  ieurc  5.6  :  Rliodamine  D,  Diffracted  Intensity  Vs  delay  ume. 


From  all  the  figures,  5.3  -  5.6,  it  can  be  seen  that  the  correlation  trace  is 
symmetric  about  the  zero  delay  position.  The  background  signal  at  r  =  is  also 
evident.  The  signal  to  noise  ratio  in  this  region  is  quite  low,  hence  the  apparent 
fluctuations. 


The  fact  of  the  correlation  trace  being  symmetric  implies  that  the  relaxation  time 
is  faster  than  the  coherence  time  oflhe  pulse,  i.e.  the  resolution  of  the  experiment.  The 
observed  signal,  therefore  gives  a  correlation  profile  |G(t)|  of  the  incident  incoherent 
light  where  G(x)  =  exp(-|i1/  Tc)  is  the  autocorrelation  function  of  the  incident  field.  In 
such  cases,  however,  it  is  still  possible  to  measure  T:  by  monitoring  both  the  dif  fracted 
beams  simultaneously  and  measuring  the  lime  lag  between  them.  For  homogeneously 
broadened  systems  the  two  traces  should  completely  overlap.  For  extremely 
inhomogenous  systems,  however,  the  two  traces  will  be  separated  in  time/16-  17 -  1<7- 
24)  This  time  corresponds  to  the  phase  relaxation  time.  Depending  on  the  degree  of  the 
inhomogeniety,  however,  this  lime  will  vary  and  is,  thus,  not  a  completely  accurate 
way  of  determining  7*;. 

Nevertheless,  the  experiment  was  performed  and  it  was  found  that  the  two 
traces  overlapped  completely.  This  could  mean  one  of  two  tilings  : 

(i)  The  materials  can  be  described  as  homogenous  systems.  From  Chapter  3,  this  is 
known  to  be  untrue. 
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(ii)  I  ! io  materials  can  ho  described  as  inhomogenotis  systems  ini l  willi  a  phase 
relaxauon  lime  .shorter  than  that  which  can  he  measured  using  the  delay  line.  The  delay 
line  moves  in  2.5 /j/h  steps  and,  since  it  is  a  double  pass  line,  this  corresponds  to  a  time 
delay  of  171s.  'litis,  then,  is  the  ultimate  resolution  of  the  experiment.  'Unis,  even  il  the 
materials  are  extremely  inhomogenously  broadened  hut  have  a  phase  relaxation  time  of 
the  order  of  20  Is,  the  two  correlation  traces  will  not  he  separated  in  time. 

Point  (ii)  is  taken  to  he  the  correct  interpretation  of  the  results  as  it  is  consistent 
with  the  findings  of  Chapter  f  It  is  also  consistent  with  theoretical  calculations  done  by 
Plyt/anis  ct  olS 

MEASUREMENT  OF  THE  MAGNITUDE  OF  THE  NONLINEAR 
SUSCEPTIBILITY 


By  measuring  the  magnitude  of  the  diffracted  signal,  when  both  the  beams  are 
time  coincident  on  the  sample,  it  is  possible  to  determine  the  absolute  value  of  the 
nonlinear  susceptibility  %'X) .  using  equation  4.17  of  chapter  4. 


(3,  _  8cV£0aVq 
1  3o)l4T(\  -  T) 


(5.19) 


The  diffraction  efficiency  is  given  by  the  ratio  of  the  diffracted  signal  to  that  of  one  of 
the  pump  beams. 


Table  5.2  displays  die  results  of  the  calculations,  including  the  results  obtained 
for  the  susceptibility  of  die  inclusions  using  the  formula  2.42. 


(5.20) 
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Sample 

Size  (A) 

Abs.  Caeff. 
[cm') 

X{i)(m2V-2) 

/, 

XY'inrV'1) 

N 

395 

221.5 

3.23x  10'19 

1.074 

3.3  lx  nr10 

M 

392 

135.6 

1.1 5x1  O'19 

1 .064 

1.98x10-'° 

() 

365 

125.3 

1.13x  10'19 

1.028 

2.34x10-'° 

S 

184 

153.8 

1.90x1  O'19 

1.330 

1.48x10-'° 

P 

52 

253.3 

7. 13x1  O'20 

1 .003 

7.86x  1()-'7 

C 

50 

50.7 

4.42x1  O'20 

1.133 

1.72x10-'° 

0 

47 

451.1 

4.95x  10'19 

1.118 

2.21x10-'° 

Table  5.2  :  Values  lor  the  macroscopic  susceptibility  and  the  susceptibility  of  die 
inclusions,  obtained  using  incoherent  light. 


X(i)  again  shows  the  familiar  square  dependence  on  the  absorption  coefficient. 
This  is  shown  graphically  in  figure  5.7.  The  values  for  £.3)  ar e  of  the  same  order  of 
magnitude  as  those  measured  by  Hache  et  al/27.  28,  29)  At  first  this  seems  surprising  in 
light  of  the  arguments  of  Chapter  3.  However,  bearing  in  mind  that  the  experiment 
measures  an  absolute  value  for  £(3),  (  as  opposed  to  a  comparative  value  )  there  is 
scope  for  error  in  the  measurement,  especially  in  the  determination  of  the  intensity.  The 
results  obtained  could  be  improved  upon  by  performing  the  experiment  relative  to  a 
well  characterised  sample,  such  as  CS2  to  allow  for  variations  in  the  beam  quality. 
This,  however,  was  not  done  as  the  results  were  only  out  by  a  factor  of  2-4,  which,  in 
nonlinear  optics  is  generally  considered  a  reasonable  error.  One  conclusion  that  can  be 
drawn  from  these  results  is  that  the  diffracted  signal  obtained  is  from  a  truly  electronic 
effect  and  not  from  a  thermal  grating.  The  logic  behind  this  conclusion  is  that  if  a 
thermal  grating  were  contributing  to  the  signal,  the  diffraction  efficiency  obtained 
would  be  orders  of  magnitude  higher. 
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figure  ,V7  :  T'*  Vs  Absorption  aK'ltlciem 

The  local  field  factor,  squared,  y, ,  at  a  wavelength  ol  522nm  is  mostly  real  and 
negative  implying  that  the  macroscopic  susceptibility  x  •  and  the  susceptibility  of  the 
inclusions  have  different  signs.  However,  in  this  experimental  configuration,  as 
in  short  pulse  DFWM,  it  is  not  possible  to  determine  the  sign  of  the  susceptibility. 

CONCLUSIONS 

The  only  conclusions  to  be  drawn  from  this  chapter  then  are  that  the  samples 
have  a  phase  relaxation  time  of  the  order  of  or  faster  than  20  (see  and  that  the 
susceptibility  of  the  inclusions  is  =  2xl()‘,"/u  \'  '  In  order  to  determine  what 
mechanisms  contribute  to  the  value  obtained  lor  x„'  it  >s  essential  to  know  its  sign  and. 

as  accurately  as  possible,  its  response  time.  The  results  obtained  from  previous 
chapters  and  this  chapter,  though  good,  are  not  conclusive  and  therefore,  it  is  necessary 
to  perform  further  experiments  to  evaluate  the  sign  and  the  temporal  response  of 

ITie  following  chapter  describes  an  experiment  that  measures  the  deflection  of  a  beam, 
due  to  the  optical  Kerr  effect,  passing  through  a  nonlinear  material  on  which  an  intense 
pump  beam  is  incident,  as  a  function  of  time  delay  between  the  probing  and  the 
pumping  beams.  The  laser  used  has  a  (see  pulse  width  and  the  temporal  resolution  is 
therefore  better  than  that  ohtained  in  Chapter  4  It  is  also  possible  to  determine,  from 
the  experiment,  the  sign  of  the  nonlinearity. 
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ULTRAFAST  BLAM  DKFLHCTION  METHOD  FOR  THH 
MFASl  JRFMFNT  OFTHF TRANSIENT  RFFRACT1 VF  INDHX  OF 

M  ATI- RIALS 


Numerous  el  feels  such  as  the  optical  Kerr  effect.' F  interferometry  and 
beam  dellection  methods,14-  can  he  utilised  lor  the  measurement  of  the  nonlinear 
index  of  refraction,//,,  in  materials.  I  he  first  two  methods  allow  time  resolution  while 
the  last,  being  a  single  beam  method,  only  permits  the  determination  of  the  magnitude 
of  iu  and  provides  no  information  on  the  dynamics  of  the  nonlinear  interaction.  If. 
however,  a  second  beam  is  used  in  an  experiment  that  measures  beam  deflection  or 
distortion,  it  is  possible  to  actually  resolve  the  interaction  in  time.  Based  on  this  idea. 
Albrecht  et  al devised  an  ultrafast  beam  deflection  technique  for  measuring  the 
induced  change  of  the  refractive  index  of  a  nonlinear  material  that  offers  simplicity  and 
sensitivity,  as  well  as  time  resolution  on  a  fs  timescale. 

The  principle  of  the  experiment  is  depicted  in  figure  6. 1.  A  weak  probe  pulse  is 
incident,  normally,  on  a  sample  placed  behind  an  aperture.  At  the  same  lime  an  intense 
pump  pulse  is  incident  on  die  same  point,  at  an  angle  a  with  respect  to  the  probe.  The 
aperture  ensures  a  definite,  nearly  rectangular  spatial  intensity  profile  and  an  accurate 
overlapping  of  both  beams  al  the  entrance  of  the  sample.  The  presence  of  the  pump 
pulse  induces  a  refractive  index  change  in  the  sample  which  consequently  results  in  the 
probe  beam  being  deducted  by  an  angle  0 .  "Hie  extent  of  the  deflection  is  a  measure  o! 
the  induced  index  change  and  its  decay,  as  a  function  of  delay  time  between  the  pump 
and  probe  beams,  gives  information  on  the  dynamics  of  the  material.  The  defection  of 
the  probe  is  measured  by  a  diode  array. 


X2 


Diode  Array 

Figure  6.1  :  Schematic  diagram  of  the  experimental  set-up  for  measuring  die 
transient  nonlinear  refractive  index  of  materials. 


The  region  of  the  sample  irradiated  by  the  pump  has  a  resultant  refractive  index 
given  by  n0  +  An,  where  «0  is  the  refractive  index  of  the  sample  in  the  absence  of  the 
pump  beam  and  An  is  the  index  change  induced  by  it.  The  change  in  refractive  index  is 
related  to  the  nonlinear  index  of  refraction,  n2  through  the  equation: 


An  =  rhI  (6.1) 

where  /  is  the  intensity  of  the  pump  beam. 

It  is  possible  to  relate  ^  to  the  deflection  of  the  probe  beam  by  using  Snells 

law. 


Figure  6.2  :  Geometrical  construction  showing  the  relationship  between  An  and 
the  deflection  angle  6. 

In  region  1,  the  probe  beam  travels  through  a  medium  of  refractive  index  given 

c 

by  nl}  +  An  with  a  velocity  V  - - .  The  total  distance  travelled  is  given  by  <1  tan  a 

n0  +  An 

.  .  ,•  d  tana, 

and  the  ume  taken  is.  therefore,  /  - ( nn  +An). 
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In  legion  2.  ilk-  probe  Ik'um  trawls  (lit i ui;_*ii  .1  nuklmni  o!  n  l  1  .icl  1  \ v  index 


will)  a  velocity  — .  I11  lime  1  u  travels  a  distance 


.1  lan  a 


•(//,,  a  .  \n  l. 


I  rom  f  inure  2.  thus,  tan  = -—tana .  When  0  —  U.  tail  (J  —*  sin  0 .  aiul  lienee 

n 


\\  e 


net: 


suit/--- — lana  10. 2) 

which  is  the  same  as  die  result  obtained  by  Albrecht  <7 

Hie  greater  the  deflection  of  the  probe  beam,  the  greater  will  lie  die  resolution  of 
the  experiment.  To  achieve  high  deflection  angles  a  is  chosen  to  be  close  to  <X6"  which 
results  in  an  interaction  length  of  =  2mm.  When  using  fs  pulses  it  is  very  important  to 
avoid  having  very  long  interaction  lengths  as  this  results  in  an  increase  in  the  pulse 
duration  due  to  group  velocity  dispersion  in  the  sample.  i;or  an  interaction  length  of 
2mrn.  group  velocity  dispersion  is  negligible.  The  diode  array  that  records  the 
deflection  of  the  probe  beam  consists  of  pixels  25(jm  wide  and  is  at  a  distance  of 
=  4l5mm  from  the  aperture.  Therefore  a  displacement  of  one  pixel  corresponds  to  a 
deflection  angle  of  6x  10'5  radians. 

The  experiment  was  performed  by  measuring  the  deflection  of  the  probe  beam 
as  a  function  of  delay  time  between  the  pump  and  probe  beams.  Hence  it  is  possible  to 
determine,  simultaneously,  both  the  magnitude  and  the  response  lime  of  the  nonlinear 
refractive  index.  It  is  also  possible  to  determine  the  sign  of  //.  by  comparing  the 
deflection  observed  for  water  with  that  observed  for  the  sample.  The  reason  that  the 
comparison  was  done  with  water  is  that  it  was  the  solvent  u.sed.  If  the  beam  through  die 
sample  is  deflected  less  dian  it  is  for  pure  water  then  n,  is  negative,  which  is  intuitively 
obvious.  This  was  the  case  for  all  the  gold  colloids  tested. 

The  laser  system  consisted  of  a  CPM  dye  laser  (  A  =  616nm,  pulse  width  = 
75fs  )  followed  by  an  Excimcr-laser  pumped  amplifier  cascade  split  by  a  beam  splitter 
into  a  strong  pump  pulse  and  a  weak  probe  pulse.  The  width  of  the  aperture  was 
1 50 flm  and  the  diameter  of  the  2  beams  was  400  )jm.  A  poiariscr  was  placed  in  from 
of  the  probe  beam,  oriented  at  45°  to  the  plane  of  the  pump  polarisation.  An  analyser 
placed  before  the  detector,  thus  made  possible  the  separate  measurement  of  the  probe 
pulse  components  parallel  and  perpendicular  to  the  pump  pulse. 
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Figure  0.3  shows  a  typical  result  obtained.  The  sample  is  a  IXnm  gold  colloid 
in  de-ionised  water.  Also  shown  on  the  same  graph  is  the  deflection  for  water  as  a 
function  of  delay  time.  The  polarisations  of  the  pump  and  probe  beams  are  parallel. 
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Figure  6.3  :  Sample  Y,  Probe  Beam  Deflection  Vs  Delay  time. 

As  is  evident  from  the  graph,  at  zero  delay,  the  deflection  for  sample  Y  is  17  pixels 
while  the  deflection  for  pure  water  is  21  pixels.  Taking  the  value  of  «,  for  water  to  be 
3.8x  10-16 cm2  /  W,  we  get,  for  the  inclusions  : 

n7  =  -7.24x  10'l7cv/t:  /  IV. 

n,  and  X  are  related  through  the  formula  :(7* 


n2  =  — 5 - 

4  n0£0c 


Using  equation  6.3,  therefore  it  is  possible  to  obtain  a  value  for  zn>  for  the 
samples  measured.  The  results  arc  tabulated  in  table  6. 1 . 
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I  able  <).l  :  lixpemncntal  Results 


Six  different  samples  were  measured,  with  eaeh  sample,  exeept  M.  being 
measured  at  two  or  more  concentrations.  As  was  observed  with  previous  experiments, 
again  there  is  no  relationship  between  and  particle  size.  However  x0)  *s  again 
found  to  be  proportional  to  the  square  of  the  absorption  coefficient.  Figure  6.4  shows 
this  result. 
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As  was  explained  earlier,  in  chapter  2.  the  nonlinearity  associated  with  the 
actual  gold  particles  is  related  to  the  effective  nonlinearity  o I  die  solution  through  : 


Where  is  the  third  order  susceptibility  of  the  gold  particles,  p  is  their  volume 
fraction  and  the  factor  ^  is  a  local  field  factor.  At  the  surface  plasmon  resonance 
frequency  it  is  greater  than  unity  and  serves  to  enhance  die  effective  nonlinearity  of  the 
solution  significantly.  Off  resonance  however,  /,  <  I  and  %t})  is,  in  fact,  diminished. 
This  can  be  seen  by  comparing  the  results  for  11  at  6l6nm  (  this  chapter  )with  those 
at  522  nm  (  Chapters  3  and  5  ). 


The  field  enhancement  factor  is  size  dependent  in  that  it  depends  on  the  width  of 
the  surface  plasmon  resonance  band,  which  should  in  turn  be  proportional  to  the  size  of 
the  inclusions.  However,  in  all  the  colloids  measured,  there  were  no  significant 
differences  in  the  widths  of  the  absorption  bands,  despite  the  disparity  in  sizes. 

/,  is  determined  using  the  equation  : 


f=  3(e,) 

‘  ei+2ed 

Where  cd  (real)  is  the  dielectric  constant  of  the  solvent  <  dielectric  matrix), 
and  £,  (imaginary)  is  the  dielectric  constant  of  the  gold  panicle. 


(6.5) 


At  A=616nm,  it  turns  out  to  be  approximately  0.7,  much  the  same  lor  all  the 
samples. 

Using  equation  6.4,  therefore,  %j3>  can  be  determined.  The  results  are  tabulated 
in  table  6.2. 
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fable  6.2  :  x\l)  values  for  the  different  dimples. 


Tlic  values  obtained  for  x\l)  are  about  an  order  of  magnitude  lower  than  those 
obtained  at  resonance.  (  A  =  522nm  ).  This  is  to  be  expected  as  the  experiment  only 
measures  the  real  part  of  the  effective  susceptibility.  In  any  case,  the  contribution  of  the 
absorptive  effect  (predominantly  imaginary)  should  be  significantly  reduced  at  this 
wavelength. 

The  decay  of  the  probe  beam  deflection  as  a  function  of  delay  lime  between  die 
pump  and  probe  beams  yields  information  on  the  dynamics  of  the  interaction.  The 
resolution  of  the  experiment  is  limited  on  the  pulse  width  of  the  laser  beam.  If  the 
response  time  of  the  nonlinearity  is  f  ast  on  the  timescale  of  the  pulse  width,  then  the 
decay  curve  is  equivalent  to  the  autocorrelation  profile  of  the  beam.  If.  on  the  other 
hand,  the  nonlinearity  responds  much  more  slowly,  then  the  deflection  decays  with  a 
characteristic  time  dependent  on  the  energy  relaxation  time,  T,  of  the  sample.  In  the 
case  of  most  of  the  samples  measured,  it  turns  out  that  7,  is  in  fact  fast  on  the  timescale 
of  the  pulse  width  which  is  75  fs.  This  is  evident  by  referring  back  to  figure  6.3.  Here, 
the  only  thing  that  we  can  do  is  put  an  upper  limit  of  75  fs  on  7,.  Two  of  the  samples 
show  a  slightly  different  behaviour.  They  are  the  two  smallest  colloids  P  and  Q  of  sizes 
of  the  order  50 A.  In  these  two  cases,  the  deflection  decays  on  a  time  scale  of  about  600 
fs.  The  decay  profile  lor  the  two  samples  is  shown  in  figure  6.5  : 
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Figure  6.5  :  Samples  I*  and  Q.  I5robe  Beam  Deflection  Vs  Delay  l  ime. 


The  reason  that  these  two  samples  behave  differently  to  the  rest  is  obviously 
related  their  small  diameter.  When  the  size  of  the  colloidal  particle  goes  below  a 
particular  threshold  level,  it  starts  to  behave  like  a  single  molecule  and  molecular  states 
come  into  being.  Molecular  states  characteristically  decay  on  a  much  slower  timescale 
than  electronic  states.  (????) 


89 


rrrrrrncrs 


<v  I  l'homa/eau.  J.,  Elchepare.  J..  Ciriilon.  G.  and  Minus.  A.  t 19X5)  ( )pi.  Leu.  10. 

225 

(>.2  Hallori,  I  ..  Okawa.  II..  Wada,  I',  and  Sisahe.  H.  ( 1902)  ( )|M.  Loll.  17.  1500 

<>.  i  lokunaga,  1-2,  Terasaki.  A..  Kohavashi.  1'.  ( 1902)  Opt.  I  .oil.  17,  1121 

0.4  Sheik  hahae.  M..  Said.  A.  A.  and  van  Stryland.  12  W.  (1989)  High 

Sensitivity,  Single  Beam  n,  Measurements.  Opt.  Lett.  14.9.55 

0.5  Sheik-bahae,  M.,  Said.  A.  A.,  Wei.  T.  H..  Hagan.  19. J.  and  van  Stryland,  E. 
W.  (1990)  Sensitive  Measurement  of  Optical  Nonlinearities  using  a  Single  Beam ,  IEEE 
J.  Quant.  Elec.  26.  760 

o.O  Albrecht,  H.  S.,  Heist,  P.,  Kleinschmidl,  J.  and  Lap.  D.V.  (1993)  Ultrafast 
Beam-Deflection  Method  and  its  Application  for  measuring  the  Transient  Refractive 
Index  of  Materials,  Appl.  Phys.  B  57,  193 

6.7  Butcher,  P.  N.  and  Cotter,  D.  (1990)  in  The  Elements  of  Nonlinear  Optics, 
Cambridge  university  Press,  306 


90 


CHAPTER  7 


CONCLUSIONS 


Hxperimental  results,  so  tar,  have  never  shown  sueh  a  si/e  dependence.  The 
reason  lieing  that  the  intrahand  contribution  (ahove)  in  actual  metal  crystallites  is  not  the 
sole  polarisation  mechanism.  Two  more  mechanisms  contribute  there  with  quite 
different  behaviour  that  dominates  that  of  the  intraband  term.  <  )ne  mechanism  is  the 
interband  term  that  arises  from  the  electronic  dipole  transitions  between  die  1  tiled  d  band 
stales  and  the  empty  confined  ones  in  the  s-p  band,  and  gives  a  contribution  that  is 
negative  imaginary  but  size  independent  as  the  d  electrons  are  unaffected  by  the 
confinement. 

The  other  mechanism  is  the  hot  electron  contribution  that  results  from  the 
modification  of  the  population  of  the  electron  states,  the  Fermi-Dirac  distribution, 
caused  by  the  elevation  of  their  temperatures  subsequent  to  the  absorption  of  photons  in 
the  resonant  process,  but  before  the  heat  is  released  to  the  lattice  of  the  crystallite;  this 
leads  to  a  contribution  to  x(3)  that  is  positive  imaginary  and  size  independent. 

From  experiments  performed  so  far,  it  seems  that  the  dominant  contribution  to 
comes  from  the  third  mechanism  and  also  that  even  the  interband  term  dominates 
the  intraband  term. 
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